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■ CDMA Patent 

c . . . 

The process at the transmitter of CDMA is roughly 

BitStream -> Frame Builder ->GoIdCode Encoding -> Modulator [Using LookUp Tablel 

D/A -> Filter! Swkched Capacitor.Clock Frequency] -> Filter! Resistors, Capasitors] «> Transmitter 

The corresponding process for CDMA receiver is 

Complex Number Seq —> CDMA Process -> Channel Estimation -> Demodulation [Using Vitcrbi] -> BitStream 
These process are in fact similar to GSM so we could build a dual mode phone. We describe the GSM processes 
The transmitter processes arc listed below 

BitStream -> Frame Builder -> Modulator [Using LookUp Table] 

D/A -> Filierf Swiiched Capacitor.Clock Frequency] Filter[ Resistors, Capasitors] -> Transmitter 

At the receiver we have the reverse process 

Complex Number Seq -> Multiplying incdmming stream by i A n -> BitSinchronisation -> Channel Estimation -> 
Demodulation [Using Vjterbi] -> BitStream 



Note 

1) Thc mode dependent process arc highlighted. 

2) The same filter can be used in both the cases if the clock frequencies are selected with care. 

3) The timing issues, and frequency dependent items like antennas have been ignored in this paper. 

4) In this paper we use many different pulses depending on the bandwidth, chip rate, required and wc trade off the non 
constant amplitude. This is the key inovation (to get higher bit rates). 

The selection of bitrate is just a mailer of engineering design involving BER, Power Amplifier efficiency and spectrum 
utilisation. We also point out that although in this paper wc have designed a modulator from first principles, we would 
in praiicc use the current Nokia Mobile Phones lookup tabic aruhitccure [ based on a pulse width of 4T ] or an 
architecure based on 6T pulse width obtaining the correct syclostatioary spectrum in each mode by the appropiaic niter 
[Analog]. The Switched Capacitor will in pratice do most of the spectral shaping common to both modes. This again is 
standard 10 design. Again the key inovation is the fact that we use one and/or two or more not n-Laurent pulses 
Co ( t ) , Ci ( ... in the construction derived through an optimisation process to construct the modulator. Others 
have used, Cq (t) only, varing at most values of B b T to vary the pulse (Linearised GMSK) , or prefiUcring to 
compensate for the reconstruction filter, wc have proposed to use two pulses (designed using an optimiser). to get lower 
amplitude variation. 

5) The key inovaiivc step is the nature of the functions used to transform the gold code into the sequence used in the 
correlator by the Function below. Il is these Functions that wc wish to protect in the patent. The transformation need Only 
occur once,and the results stored, as the mobile can pre calculate the sequence d t for i = 0,1 , .. N-l given the code 
( c o. cj c N _i} when the code is asigned. 



Transformation 1 (to detect +1 symbol, C 0 equivalent 
pulse ) 
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yi = fori = 0. I. 2 N — I 

Given code ( q>, c\ , cn_i 1 where N is the number oF elements in the sequence. 

3j = 1 if Ci = I ; for i = O, ... N - 1 
^ = -1 if d =0; fori = 0, ... N — I 
^ j — ao ; 

bj = bi_, + aj Fori = 1, 2 N — 1 ; 

dj - y» for i = 0, 1 , 2 N — 1 and 

• = V(-1) 

There is another transformation that can also be used 

Transformation 1 b (to delect + I symbol, Co equivalent pulse ) using the same notation for di 
di = r b > fori = 0, 1. 2,".... N — 1 and 

Transformation 2 (10 detect -I symbol.Co equivalent 
pulse ) 

yi = fori = O, 1, 2, .... N - 1 

Given code ( c 0 , Ci , c N .j } where N is the number of elements in the Sequence . 

aj - -1 if a = 1; fori 0, ...N - 1 
aj = 1 if ci =0; fori = 0, ...N— 1 
bo = ao ; 

bs = bi_i aj for i = l, 2 N - 1 ; 

di = yi fori ^ 0, 1, 2 N — 1 and 

' = V(~n 

There is another trans formation that can also be used 
Transformation 2 b (to detect — 1 symbol, Co equivalent pulse ) 
d; = i" b - fori = 0. 1. 2 N - I 



Transformation 3 (to detect +1 symbolC^ equivalent pulse ) 
yi = (-I)* fori = 0, 1. 2 N « 1 

Given code { Cq , c j , e N _j ) where N is the number of elements in the sequence . 

a; = 1 if cj = 1; fori = 0, ... N - 1 
ai = -1 if Cj = 0; fori = 0, ... N - 1 
bo = ao — /-t- ajsi ; 

bi = bi_i + a» - ftj.i fori = 1, 2, .... N - 1 ; 
di = yj i b > fori = 0, 1, 2, .... N — 1 and 
i = 

There is another transformation that can also he used 
Transformation 3 b (to detect 1 symbol, Cj equivalent puUe ) 
dj = r*> fori = 0, 1, 2 N- 1 and 

Transformation 4 (to delect — I symbolCi equivalent 
pulse ) 

yi = (-1) 8 fori = O. 1,2 N — I 

Givencode ( c 0l c h c N -j } where N is the number of elements in the sequence . 

fli = -lif ^ s 1; fori = 0, ... N — I 
a; = 1 if cj = 0; fori = 0. ...N- 1 
bo = ao — /+ aN ; 

bi = bi_ t + a 4 - ai_! fori = 1,2 N — 1 ; 

di = yi t b > fori = 0, 1, 2 N- I and 

, = V(-l) 

There is another transformation that can also be used 
Transformation 4 b (to detect - 1 symbol, d equivalent pulse ) 
di » for i = O, I. 2 N - 1 and 

We also claim any 

a)cyclic rotation of the sequences resulting from Transformations 1.2.3.4 and lb.2b,3b.4b 
^multiplication of the elements of the sequence obtained from Trans formations 1,2.3.4 and lb.2b 3b 4b 
constant (real or complex) 
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are also claimed to be transformations that we wish to patent. 

We also claim any combination of cyclic rotation followed by linear scaling of the result is patented. 
-"Hio subroutine CodingTransformation implements the transformations 1.2,3,4 



We also stress that while the bitstream is differentially encoded and the substitution used is {0 -» -l } the 
substitution elsewhere used for the gold codes is (0 -* — ], 1 -> 1 ] 

6) The demodulation in the CDMA mode uses the 4 separately transformed codes( Transformation 1,2 3 4) or 
(Transformation 1 ,2). Also we point out that in this scheme we have two adjacent chip positions that generate | arg 
values of correlation. With thus using 2 pulses and 8 complex correlations, we can obtain, should we desire the best 
possible demodulation. Thus, most of the work is done by the gold code, during correlation. 

7) There is nothing special about the Gold Code used in the example. Any gold like code will do as long as it posscscs 
the usual properties of a gold code. 

8) In case we consider a BT = 0.15 etc GSM modulation scheme, ( i.e. twice the speed) we can keep the 

Example of the output using different encoders which encode transformations 
to correlate with the input of a bipolar bit stream of the first training sequence of 
GSM system where we have used the rule {0 -1} and not differentially 
encoded the . 

T ^f ri ^V^ ining ? cqucnc * of GSM s y stem whcrc wc ^vc used the rule [0 -> -1 ) and not differentially encoded it is 
eallcs GSMseq ={-1 .-1,1,-1,-1,1,-1,1,1,1 -1-1-1 -1.1,-1,-1-1.1-1-1 ,1 ,-1,1.1.1) 

ModOutputGSM = Modulator [TO [ GSMaaq, NumberOfCurves 2, 

ModulAfcingrPuaae Opt Pulse Sea led, Samplinglnte^l -* T / 4 ] ; 

We call the Modulator output ModOutputGSM. PrimitiveCDMARcccivcr encodes GSMseq according to 
transformation rule 1 -it, 

ToznGSMl = PrimltiveCDMARecaiver[ModQUCp\ltaSM, 1, 4] j 

LidtPlot^wnGSHl/ZAba, Plot Joined -> *rue, FlotRange -> All] 




Graphics 



We have plotted die autocorrelation. The next graph shows the autocorrelation when GSMseq is encoded with 
transformation ib. It seems the same as with transformation lb 



PrijnitiveCDKAIlo C aiv a r2 [WodOutputGSM, GSMaeq, X, 4]; 
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LiatPlot[% // Aba. PlorCToined -> True, flotRange -> All] 




- Graphics - 

PriminveCDMARcccivcrGSM2Pulse encodes GSMseq with Transform 3 

Pzrlxo±t:±veCX}HJUlecei.vQxrGSH2Pulfle[KodOutpvifc6SM. GSMseq, 1, 4]; 
tiptPlot [«6 // Abs ( PlotJoined -> True, PlotHanga -> All] 




- Graphics - 

PrimitiveCDMAReccivcrGSM2PulseEfficieni encodes GSMseq wkh Tmnsform3b. This seems io be diffrcnt 

Pi:imlCiveCDMAR©ceiv©rGSM2P\ilaGE££ic:iG»t[ModCait:puCGSM # GSM&aQ;, l r 4]; 
tiistPlot [9* // Abe, PlotJoined -> Trua, Plot Rang© -> All] 




- Graphics - 
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■ Dual Mode CDMA and GSM Operation 

in this work we demonstrate how lo build a dual mode phone using as much as possible the same hardware for both the 
systems. 

Given the following process for GSM Transmitter 

BhStream Frame Builder -> Modulator [Using LookUp Tabic] 

D/A -> Filter[ Switched Capacitor.CIock Frequency] -> Filter[ Rcsistors,Capasitors] -> Transmitter 



At the receiver we have the reverse process 

Complex Number Scq — > Channel Estimation -> Demodulation [Using Vjrerbi] — > BiiSircam 
The corresponding process at the CDMA transmitter is roughly 

BitStream -> Frame Builder ->GoldCodc Encoding -> Modulator [Using LookUp Table] 

D/A -> Filter[ Switched Capacitor.Clock Frequency] -> Filler! Resistors,Capasitors] -> Transmitter 

The corresponding process for CDMA receiver is 

Complex Number Scq — > CDMA Process -> Channel Estimation — > Demodulation [Using Binary Decision] — > 
BitStream 

Note 

1 )The mode dependent process arc highlighted. 

2) The same filter can be used in both the cases if the clock frequencies are selected with care. 

3) The timing issues, and frequency dependent items like antennas have been ignored in this paper. 

4) In this paper we use many different pulses depending on the bandwidth, chip rate, required and we trade off the non 
constant amplitude. This is the key inovation (to get higher bit rates). 

The selection of bitrate is just a matter of engineering design involving BER, Power Amplifier efficiency and spectrum 
utilisation. We also point out that although in this paper wc have designed a modulator from first principles, we would 
in pratice use the current Nokia Mobile Phones lookup table architccure [ based on a pulse width of 4T ] or an 
architccurc based on 6T pulse width obtaining the correct syclostatioary spectrum in each mode by the appropiate filter 
[Analog], The Switched Capacitor will in pratice do most of the spectral shaping common to both modes. This again is 
standard to design. Again the key inovation is the fact that wc use one and/or two or more not n-Laurent pulses 
Co ( t ) , Ci ( c) ... in the construction derived through an optimisation process to construct the modulator. Others 
have used, C 0 < t ) only, varing at most values of B b T to vary the pulse (Linearised GMSK) > or prcfiltcring to 
compensate for the reconstruction filter, we have proposed to use two pulses (designed using an upiimiser), to get lower 
amplitude variation. 

5) The key inovative step is the nature or the functions used to transform the gold code into the sequence used in the 
correlator by the function below. It is these functions that we wish to protect in the patent. The transformation need only 
occur once.and the results stored, as the mobile can pre calculate the sequence di for i = 0.1. .. N-l given the code 

I cq . Cj cn_i } when the code is asigned. 



Transformation 1 (to detect +1 symbol. C 0 equivalent pulse ) 
yj - (-1)' fori 0. 1, 2 N-l 

Given code [ c 0 , c t . ..... c N _j } where N is the number of elements in the sequence . 

a; = 1 if Ci - 1; fori = 0, ... N — 1 
ai = - I ir cj = O; for i = O, ... N — 1 
bo = a© ; 
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bj = + ai fori = 1,2 N - 1 ; 

di s y L fori = 0, 1, 2, .... N — I and 

X~ V(-i) 

ihere is another trans formal ion that can also be used although it does not perform as well 
Transformation 1 b (to detect + 1 symbol, Co equivalent pulse ) 

dj = i _b ' fori - 0, 1, 2 N - 1 and 

Transformation 2 (lo detect —1 symhol,C 0 equivalent pulse ) 
yi = (- 1) 1 fori = 0, 1, 2, .... N — 1 

Given code { c 0 , Ci , Cn -i ] where N is the number of elements in the sequence . 

a; = — 1 if Ci = 1; fori - 0, ...N- 1 
aj = 1 if Ci = 0; fori = Q t ... N - 1 
bo = ao 

bj = bj-i + a; fori = 1, 2 N - 1 ; 

di = y; fori = 0, 1 , 2, — N — 1 and 

i = Vc-n 

There is another transformation that can also be used although it does not perform as well 

Transformation 2 b (to detect - 1 Symbol, Co equivalent pulse ) 

di = i-* fori = 0,1,2 N - 1 and 

Trans rormau on 3 (to detect -H symbolCi equivalent pulse ) 

yi » (-1)* Tori « 0, 1, 2 N - 1 

Given code { cq , Cj c N _| } where N is the number of elements in the sequence . 

a 4 = 1 if c ; = 1; fori = 0, ... N - 1 
a s = -1 if Ci = 0; fori - 0, ... N - 1 
bo = no — /•*- 3n ; 

bi = bi_i + ftj - aj_! fori = 1,2 N — 1 ; 

di « yi i bl fori = 0, 1, 2 N - 1 and 

i = V(-i) 

There is another transformation that can also be used although it docs not perform as well 
Transformation 3 b (to detect 1 symbol, C, equivalent pulse ) 

di - fori = 0, 1. 2 N - 1 and 

Transformation 4 (to detect -1 symbolCi equivalent pulse ) 

yi = (-1V fori = 0, I, 2 N — l 

Given code [ Co, c, cn-j I where N is the number of elements in the sequence. 

a» - - J if q = 1; Tori = 0, ...N - 1 
a ; = 1 if cj = O; fori = 0, ... N — 1 
b 0 = ao - /+ a N ; 

bj = bj-i -»- ai - aj-i for i = J, 2, .... N — 1 ; 

di = y» fori = 0. 1, 2 N — 1 and 

i = V(-l) 
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There is another transformation that can also be used although it does not perform as well 
-Transformation 4 h (to detect - 1 symbol, C t equivalent pulse ) 

si r bi for i = 0. 1,2 N — I and 

We also claim any 

a)cyc|ic rotation of the sequences resulting from Transformations 1.2.3.4 and lb.2b,3b,4b 
^multiplication of the elements of the sequence obtained from Transformations 1.2,3,4 and lb,2b.3b,4b by a 
constant (real or complex) 

arc also claimed to be transformations that we wish to patent. 

We also claim any combination of cyclic rotation followed by linear scaling of the result is patented. 
The subroutine CodingTransformaiion implements the transformations 1,2,3,4 



We also stress that while the bitstream is differentially encoded and the substitution used is (0 ->1 , 1 -> ) the 
substitution elsewhere used for the gold codes is (0 -* -1, t -» I } 

6) The demodulation in the CDMA mode uses the 4 separately transformed codes( Transformation 1.2.3.4) or 
(Transformation 1,2). With thus using 2 pulses and 4 at most complex correlations, we can demodulate the bitstream. 
Thus, most of the work is done by the gold code, during correlation. 

7) There is nothing special about the Gold Code used in the example. Any gold like code will do as lone as it posseses 
the usual properties of a gold code. 

8) In case wc consider a BT = 0.15 etc GSM modulation scheme, ( i.e. twice the speed) we can keep the 

Example of the output using PrimitiveCDMAReceiverGSM2Pulse which encodes transformation 2 to correlate with die 
input of a bipolar bit stream of the first training sequence of GSM system where wc have used rhe rule (0 -> -1) and 
not differentially encoded the . 

The first training sequence of GSM system where we have used the rule f 0 -> -1 ] and not differentially encoded it is 
calles GSMseq =[-1,-1,1 -1-1. i-ijjj .^1,-1,-1,-1.1.-1.-1 -l _i, 1,-1.1,1,1 ) 

ModOutputQSM = Modulator [L] [ GsHseq, Numberofcurvefl ^ 2. 

ModulatingPuleo QptPulaeScai«d r S^nnp ling Interval T / 4 ] j 

We call the Modulator output ModOutputGSM. PrimiuveCDMARcccivcr encodes GSMseq according to 
transformation rule I 

TomGfiwl = PrimitivaCDMftReceiver [ModOutputQSM, GSMdeax, l, 4] ; 
LiatPlot[TomGSMl // AJbs, PlotJolned -> True, PlotR^uige -> AX1] 




- Graphics 
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Wc have plotted the autocorrelation. The next graph shows the autocorrelation when GSMseq is encoded with 
transformation lb. It seems the same as with transformation lb 

t»rimitiV6CDMAReceiver2 [ModOutputGSM, GSMseq, 1, 4] ; 

Li at Plot [% // Abs, PlottToin^a -> True, FlofcRange -> All] 




- Graphics - 

PrimitiveCDMARcceiverGSM2Pulsc encodes GSMseq with transform 3 

Pri^tiv«CI3MARGceiv©rGSM2Piilae[ModOutputaOM # GSMseq, 1, 4]; 
LifltPlot[%/yAbfl. PlotJoinod -> True, PlotRango ^> All] 




- Graphics - 

PriinitiveCDMARccciverGSM2PulscEfficient encodes GSMseq wiih iransform3b 



6rimitiVBCI3MAReceiv6i:GSM3l»ulaaBf£i 



ciTOtLModOUfcpufcGSM, GSMseq, 1, A]; 
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ListPlot[^ // AbSy Plot Joined -> Txruo, FlotRcwge -> All] 



c 



10 
8 
6 
4 






( 10 


nJ 15 


V20 


V 25 



- Graphics - 
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■ GSM with BT = 0.15 

We propose u«ein« the same frame structure without differentially encoding the bits. We make the transformer, {0 -> 
Wc describe the current GSM processes 

'JZ^^S^^ "*crontialI y encoding the data but keeping the same frame suture. Wc use the current or 



u !2 in ^ir c t™!^r ! VC ,S ^ ° f ^ fUnC " 0nS USCd t0 fining sequences into the sequence 

Z^fl C °/ Telatl ° n l ° the .mpulse resposc in the Viterbi. They are listed below. It Is these functions thaiS 

ta^i^? ^orT-^ lT T r ?™* l0 l nCCd ° n,y ° CCUr ° nCe ' nnd thB results — d ' - mobile canVf i^l^ 
inc sequence di for 1=0,1,.. N-l mven the iraininfiscauence It, n. ^. .1 u kt _ ^ F 



Transformation 1 (to detect +1 symbol, C 0 equivalent 
pulse ) 

yi = (-I) 4 Fori = 0, J, 2 N - I 

Given code ( c 0 , c , , c n-i } where N is the number of elements in the sequence . 

&i « 1 if Ci = 1 ; for i = 0. ... N — 1 
a; — — 1 if Cj = 0; fori = 0, ... N — 1 
bo = ao ; 

- b^i + aj fori = 1,2 N — 1 ; 

di = yj i K for i = 0, 1 , 2 N - 1 and 

i = V(-l) 

There is another transformation that can also be used 

Transformation lb (to detect - 1 symbol, Q, equivalent pulse ) using the same notation for * t 
di = i * fori = 0, 1, 2 N- 1 and * 

Transformation 2 (to detect -I symbol,c 0 equivalent 
pulse ) 

y; ^ C-IV fori = 0, i, 2 N — l 

Given code ( c 0 . c, c n-i I where N is the number of elements in the sequence 

ai = - 1 if cj = 1 ; for i = 0, ... N - 1 
fli = 1 if = 0; fori - 0, ... N — 1 
bo = 3q ; 

bi - bj-| + a; fori = 1, 2, .... N - 1 ; 

= yi i b > fori = 0, 1. 2 N - l and 

i = V(-l) 

There is another transformation that can also be used 
Transformation 2 b (to detect - 1 symbol. Q> equivalent pulse ) 
dj = i * for i - 0, 1 , 2 N - 1 



Transformation 3 (to detect +1 symbol^ equivalent pulse ) 
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y, = (~\? fori = 0. I, 2 N — 1 

Given code (c 0 . c, , c N ,, J where N is the number of elements in the sequence 

a; = 1 if Cj = I ; for i — 0, ... N — 1 
^ »i = - 1 if Ci =s 0; for i = 0, ... N - I 
V o = ao - /+ a N ; 

bj = bi_, + aj - a ; _, fori = 1. 2 N - 1 ; 

d. = y* i bi fori = 0, I. 2 N — J and 

' = V(-l) 

There is another transformation that can also be used 
Transformation 3 b(Lo detect 1 symbol, Ci equivalent pulse ) 
d; = i fori = 0, 1. 2, .... N — | and 

Transformation 4 (to detect -1 symbolC x equivalent 
pulse ) 

yi = fori = 0, 1, 2 N — 1 

Given code { c 0 , c, c N ) where N is the number of elements in the sequence 

&i = -1 if C; = 1; fori = 0, ... N — 1 H 
a j — 1 if = 0; Tori = 0, ...N - I 
bo = tu> — /•*- a N ; 

b ; = bi_, + ai - ai ,, fori = I, 2, .... N — 1 ; 
* = ^ i b > fori = 0, 1, 2, .... N - 1 and 
i^V(-l) 

There is another transformation thai can also be used 
Ti^nsforxnation 4 b (to detect - 1 symbol. C, equivalent pulse ) 
o\ = i * fort r* O, 1. 2 N— 1 and 

We also claim any 

a)cyclic rotation of the sequences resulting from Transformations 1,2,3,4 and lb 2h 3b 4b 
constant^ Transformations' 1.2.3,4 and 1 b,2b.3b.4b by a 

arc also claimed to be transformations that we wish to patent 

We also claim any combination of cyclic rotation followed by linear scaling of the result is patented. 
The subroutine CodingTransformation implements the transformations 1,2,3,4 

^ZS^^J^ thC r H S ,- ° f ^ bitStream iS d * f *™*"y encoded and the substitution used is { 0 I ^ _ U 
the substitution elsewhere used lor the training sequences is {0 -» -1, 1 - ] } ' 

6)Thc demodulation in the CDMA mode uses the 4 senarafr»fv .r.n.F^^ ^ ^ / -r 

^TransfnrmAtmn i o\ \xr -.w u v „ L * se P arate *y transformed codes( Transformation 1,2.3 4} or 

2^2 ZSSSS^SXZ™* Cod * used in the examplc ^ *** «- - - >o« 8 as it possoS e S 

8) In case we consider a BT = 0. 15 e,c GSM moduladon scheme. < i.e. twice rhe speed) we can keep the 
»o fl u lacla9Pul8a _ QptPul- e3ca ied: a «nplinS™ a r^ T / 4 ] j 
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We call the Modulator output ModOutputGSM. PrimiiiveCDMARecei 
transformation rule 1 



ver encodes GSMseq according to 



Q Ton^SMl = PrimitlveCDMAReceiver [KodOutputQHM, GSMaeq, 1, 4] ? 

ListPlottToaGSMW/Abs, Plot Joined -> T «i e , PlorR*^ -> A i X] 




10 U 
- Graphics - 

^ e „! , f«™ P ! 0 " C, ?u h ? nulocon ; e,alion - The "«t graph shows the autocorrelation when GSMseq is encoded with 
transformation 1 b. Ii seems the same as with transformation lb enuoacd witu 

er-Jjnlti^CM«ARecei.»er2[ModQ«tput:OSM, GSMseq, i, 4 jj 
t±sfcPlot[%// Abfl, SlotJMnBd -> Trus, PlotRan^e -> All] 



f 2 
ao 

A / 8 

V 6 






10 


l5V-> 


f 20 


25 



- Graphics - 

PrimiiiveCDMARcceiverGSM2Pulse encodes GSMseq with Transform 3 

PxiJiativoCXIMARaceivexoeMaPuaoeCMoaOutinjiteSM. GSH H e q< 1, 4]; 
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t*iatPlot [*6 // JOdh, PlotJolxied -> True, Plot Range -> All] 



c 




- Graphics - 

PriTr.itiveCDMARcc E iverGSM2PulscEl*ncientencodes GSMseq wiLh Transform3b. This seems to be diffrcnt 

PriadtiVACDKARaceiverGEMJlKilaeE^ficientCltodOllt^ GSMaeq, 1, A]; 

ListPlofc[* // Abe, PlotJolned -> True, PlotRanga -> A il] 
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Needs [ M Ri.ngFuxLce±oziaD:l.££&xicoaea* - ] 
N&naa [ " RingPunc t ioaeDi f ffincoded * * - ] 

Ringpoweir. Sequen ce Generat„. Specif fedGold^mU. 
Needs [ "LanrentFimctiorifl * » j 

RuleDalayed: :rhS : Pattern t_ appears On Che right-hand aide of rule 

PhaseAnglef^] [c_J (Pha S eAngl e[ L] r t _j = Module [ {x^ *2. *1 . x< ^ x6) . 

Needs [ -LaurenttJOtationTest * ■ ] 

NeSd3: :nOCOnt : Context LaurencNotAtionTesc' was not: creatad when Needs was ev.lua.ed. 

Information on the Functions used can be obtained using help. 
Names [ "LaurentFonctiaae * * ] 

(AXN, AlphaKX, ANKInicialStateSetUp, BT, FiltPulfie, h, hFi Itemed InitialStar^ T 

?^?f?£?^ V ^', PhaSeAn ? le ' PhaseA ^5XeFasc, Receiver. Receiver Proper S A * C ° r ' 
Samplmgmterval, scartingQuadr^int, SyncSample. T. C, *, ^> ' 



* := 

912500 
:= 0-3 



Modulationlndex 2 = 

2 

<< Modul at orDat a . m; 



<< OptiinalPulfleShapefl .m; 
Plot [Opt»uXee[t.J [a) [tj. {fc, 0 




Graphics 



Tatble [ 
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2 



c 



Plot[Opcmee[L] [1] [t] , {fc, 0.6}] 
0.2r 



0.15 




812500 

The unit of time is T = 1 for OpiPulse.We scale the Pulses 10 T = 



for the unit of time. 



OptPulfleScale<a[8] [0] : = Opt^ulse [1.] [O] [tVT] 

QptPulflB3calad[B] [1] [t_J r = OptPulfi«[L] [1] It / T] 

Plot [OptPulpoS C aled[L] [D] [fc], {fc, O , BT>] 




SxlO"* 3 • Q0DaD:LOO00c X5000flD200OOQ5OQ003 
Graphic si - 



{1. 1, -1. -1. -1. 1. 1. -1, 1, -1. l. „1, 
-1. 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, _1, 
1. -1. 1, -1. 1, -1. 1, 1, -1, -1 , 1, -1, 
1, 1, 1, -1. 1, -1, -1, 1, 1, 1. -1, -1/ i. 



1. -1. 
1, -1. 
1. -1. 



-1. -1. -1, 1, 1. 1 # -1, -1. -1, 1, 1. 1, 

1. -1, 1. -1, -1, 1, 1, 1, 1,-1, 

1, 1, 1, 1, -1, -1, -1, 1, -1, 

1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, -1} 
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Plot [F±3.tPuls«[Ii] [0] [t] r {fc, O, Btt)] 



c 




5x10"®- 000© 100 00a 500 0© 20000(250 00 03 
- Graphics - 



m The Gold Sequence Set 

We generate the sequences using the method specified by Scrdar Boztas and P Vijay Kumar in Ref [1], The numbering 
of the sequences is the one used in the paper. We generate a small subset of the sequences. There are 2 10 + I sequences 
with the quaternary polynomial used. Given any binary primitive polynomial, we can generate the corresponding 
quaternary polynomial. 



Golds enlist = 

Specif iedGoldSoquencafi [ {1, 3, 2, X, O, 3, 0, 0, 1)][{1, 2, 3. 
The last sequence of scqlist has nice autocorelation properties. 
The last sequence is in fact a rn-scquence of length 1023 bits. 



t - a. Last} J i 



Goldseallsfc // Last 



{1, 0. 
0, 1, 

1- 1* 

0, 1, 



0. 

1. 
1. 



1, 
1. 
1, 



1. 1 

0, 1 

1, o. 



0, 0, 

0, 1, 

1, 1. 

0, 



Or 0. 

0, 1. 



0» 1. 
1. o. 0, 
1, 1. 0. 1. 1. 



0, 1, 1, 1 

1. 1. o 
1, o, o 



0, 

0, 0, 

1. o; 



0, 1. 

0. 1. 

1. 1. 



1, 1- 



1, 1- 1, 0, 



1, 0. 0. O, 1, 



1, 


1. 


0. 


0. 




0, 


o, 


o. 


1, 


Q, 


o. 


1, 


1, 


0, 


1, 


0, 


0 , 


o, 


1, 


0. 


1, 


1, 


1, 


1 , 


o. 


0, 


1, 




0, 


0, 


1, 


1. 


1, 


0, 


o. 


1. 


1. 


1. 


0, 


1, 


0. 


1, 


1, 


o* 


1. 


1, 


0, 


1. 


0, 




1, 


0, 


0. 


0. 


1. 


1, 


0. 


1, 


1, 






1, 


1, 


1, 


0, 






1, 


1. 


1. 


0, 


1, 


1. 


1, 


1. 


1, 


0. 


1, 


1. 


1, 


o. 






o* 


1. 


0. 


1, 


0, 


0. 


0. 


1- 


0. 


1. 


0. 




1, 


1, 


1. 






o, 




1, 


0, 


1, 


0, 


0, 


0. 


0. 


1, 


0. 


0. 


0. 


o. 


1. 






o, 




0, 


1. 


0, 


0, 


1, 


0, 


1, 




0. 


1, 


1, 


o. 


o r 






1, 




1, 


1, 


0. 


0, 


0, 


1* 


1, 


1, 


1, 


1. 


1, 


1, 


1, 




1. 



AutocorralationSeguencefGoldfleqliqt / / I*ast] 

{255, -1, -1, -l. -i. .1, _i, -1. -1. -i, _ 1# _i, _ 1# _! 

-1. -1. -1, - 1, -1, -1, -1, -l, -1, -l, -1, -1, -1, _ 1( 

"J- "1- -1* -1- -1, -1, -1, -1, -1, -1, -1, -l, -1, 

-J' -}> "J- "1* -1* -1, -1, -l, -1, -1. -1, -1. -1, -1, 

"J' ^' "3" ~ 1 - - 1 - - 1 ' - 1 ' 

"J' - 1 - - 1 - - 1 . -i. -i. -i, -i. 

-1, -1* -1, -1. -1, -1, -1, -1, -1, -1, -l, -l, 

- 1 - - 1 - -i. -i. -1. -1. -1, -1. -1. -1. -1, -X, -1, -1, 

"i # 'I' ~J' ~ X ' _1 ' _1 ' - 1 ' - 1 ' 

-1. -1. -1, -1. -1, -l. -l, -l, -1, -1, _x. -1, -i, -i, 

-1, -1. -1, -1, -X. -X. -1, -1, _l ( -1, _x, _x, -i, -x, 

"I* "I* -1. -1. -1, -1. -1, -1, -1. -1. -X, -X, -1, 



-1 
-1. 

-1, 
- 1 , 
-1. 
-1, 
-1, 
-X, 

-1, 
-1, 
-1, 
- 1 , 



-X 

-1. 
-1, 
-x, 
-1, 

-It 
-1. 

-1. 
-1. 
-1. 

-1, 



-1 
-1. 
-1. 

-1, 

-X, 

-1. 
-1. 
-1. 
-lr 
-1. 
-1, 



-1. -1. 

-1, -1, 



"lr 
-lr 
-X, 
-1, 
-1. 
-lr 
-1. 
-1, 
"I, 
"lr 



-1 , 
"lr 

-X, 
-1. 

-1* 
"I, 
"I, 
-1, 
"lr 
-1* 



"lr ' 

-lr - 

-1, - 

"lr ' 

"lr " 

■1. - 

1, - 



lr "lr 
1- 



X, 

lr 
lr 
lr 
1, 
1, 



-1. -x. 



■1 r 

"lr 
■1, 

-1, 
-lr 

-1, 



1} 



The third sequence in the list has the following autocorrelation 
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Autocar r«lat±onSeciii«nc© [Goiaoefjliet / / # [ [3] ]6] 

{255, -1. -1, -1. -1. -1. -1, -17, -1. -1 ( -1, -17. -1. -17, -17, -17, -1, -17, -1, -1. 
-1, -1, 15, -1, -1. 15, 15, -17, -17, -1. -17, 15, -1, -1, -17, 15, -1, 15, -1, 15, -1, 
15, -1, -1, 15, -17, -1, 15, -1, -1, -17, -17, -17, -1, -17, 15, 15, 15, -1, -1, 15, 
-17, 15, -1, -1, -1, -1, 15, 15, -1, 15, -1, -1, -17, -17, 15, -1. -1, 15. 15, -1, -1, 
-17. -1, -1, -17. -1, -1, 15, -I, 15, -17. -1, -1, -17. -1. -1, 15, -1, -17, 15, -1, 15, 
-1. 15, 15, -1. 15, -17, 15. 15, -1, 15, -1, 15, -17, -1. -1, -1. 15, -17. -17, 15, -1. 
-17, -1, -1. -1, -1. -1, -1, -17, -1, 15. -17. -17, 15, -1, -1, -1, -17, IS, -1, 15, -1, 
15. 15, -17, 15. -1, 15. 15. -I, 15, -1. 15. -17. -1, 15, -1, -1, -17, -1, -1, -17, 
15. -1, 15, -1, -1, -17. -I. -1, -17, - 1 . -1, 15, 15, -1. -1, 15, -17. -17, 
-1. -1, 15, -1, 15. 15, -1, -1. -1. -1, IS. -17. 15, -1, -1, 15, 15. 15, -17. 
-1, -17, -17, -17. -1, -1, 15, -1. -17, 15, -1, -1, 15, -1, 15, -1, 15. -1, 
15, -17. -1, -1. 15, -17. -1, -17. -17, 15, 15, -1. -1, 15, -1. -1. -1, -1, 
-17, -1. -17, -17, -17, -1, -17, -1, -1, -1, -17. -1, -1, -1, -1, -1, -1} 



We generate the output of the modulator 



HoctOutput: n Modulator[L] [ (Goiaeeqliflt // Laet) / . {0 -> -1) „ Number Of Curves -* 2. 
MoaulAfclngpulse -> Filtpuloe, Samplinglnterval -* t / 4] s 



V/e check the output 



ListPlot [ (Ra tModOutpvLfc] , Ln [HodOutpub] > // Transpopa, Plot Joined -> True, 
ASpectRatio -> 1] 




- Graphics - 



Now we try to test the modulated sequence using the 



The number of samples per chip is equal to 




1 
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c 



(RacaivarfL) [WodOutput, StarfcingQuadrant -> o, 

SamplingintorvaX -» T/4, HodulatingPulsft -» FilfcPulsel //Drop[« f 41&5 - 
C (Goldaeqliet // L^Bt // Drop[#. -4]&) /. {0 -> ^JL)> 



0, Q , 
0 ( 0, 
o# 0. 
o. 0, 

0. 0. 
o, 0, 



0. 0. 



0, 

0, 
0, 
0, 



0, 0, 



0. 
0, 



0, 0, 
0, 0, 



0, 

0, 



o* o, 
a, o. 



0. 0, 

o, o, 



0, 0, 
0, 0. 



0, o, 
0. 0. 



0 
0, 

0, 



0. 
0 



0, 0. 
0. 0. 



0. 0, 
0. 0. 
0, 0, 
0, 0. 

a, o, 
0, o f 

0, 0, 
0, 0, 
0, 0, 



0, 0, 

o. a , 

0, 0, 
0, 0. 
0. 0, 
0. 0, 
0. 0. 
0. 0, 
0. 0, 



0, o, 
0. 0, 



0, 

o ( 

0. 



0, 0 
0 
0 



0, 0, 
0 . 0. 



0 

a, o, o 

0, 0, 0 
0, 0 

a, o, o 

0, 0. 0 



o, o, a, o, o, o, o, o, 

O, 0, D, 0, 0, 0, o, 0, 



0, O, 0. 0, 0, 0, 0, 

o. o, o, o, a, o, o, 

0. 0. 0. 0. 



0, 0, 0, 0, o, 

0, 0, 0, 0, 0, 

0, 0, 0, 0, o, 

0, 0. 0. o, 0, 

0 , 0. 

0. 0. 0, 0, 0, 



o, 
o, 

0, 
0. 



0, 0. 
0. 0, 

0, 0, 0, 0. 0 , 0) 



Wc have successfully demodulated the biistream 



ModOutputOpt = Hadulator[L] [ (Geldseqllsc // Last) /. {0 -> -i>, 

Numbarofcurved -> 2, ModuiacingPuls© -+ OptFulseScaled, sampllnsrintenrai -+ T/4]! 

ListPlot [ (R« [Modoubputopt] , Im[MDdOutputC>pt]} / / -Pjranspoe©, HotJoined -> t™, 
AspoctRatio -> 1] * 




Graphics - 



i /« 3 ^ ,lil i» Illt f rv ^l ~> t/4. Modulatlngpulsa- optFulse sealed] //Dropr# r 4]*) 



{0. 0, 
0, 0, 

o, o, 



0 

0, 0 
0, 0 
0. 0 
0, 0 



0. 



0. 0, 
0. 0, 



0, 
0, 



0, 0, 
0, 0, 
0, 0, 



0, 0, 
0, 0, 



0, 
0, 



0, 0, 
0, 0, 



0, 0, 0, 0, 



0, o, 
o, o, 
o, 0, 
0, o, 



0 
0 
0, 



0, 0, 
0, 0, 

0, 0, 



o, 0, 
0. 0, 
0, 0, 
0, 0, 
0. 



0. 0. 



0. 
0, 



0. 
0, 
0, 
0, 



0, 0, 
0, 0, 



0. 0. 

0, 0, 

0, 0, 

0, 0, 

0. 0, 

0. 0, 

0. 0. 



0. 
0. 
0. 
0. 
0. 
0, 
0, 
0. 
0. 



0, 0. 0, o, 
0. 0, 0, 0. 



0 

a 

0. 0 
0, 0 



o, a, o, o 

0. 0, 0, 0 

o, o, a 

0, 0, 0 
0, 0. 0 



0, 0, 



0,0.0, 0, 0, 0, 

0, 0, 0, o, 0, 

0, 0, 0, 0, 

0. 0, 0, 0, 

0. 0. 0, o, 

0, 0, o, 0, 

0, 0, 0, 0, 

o. o, o. 0, 



o, o, o, o, a, o, o, o) 



■ CDMA OPERATION 

As an example given a symbol stream e.g. [1,1,-1,-1, 
demonstrates the encoding process. 



} consisting of -I .and 1, ihe following function 
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Moaule[{xl, x2, x3) # 

Xl e GoldSflQ / - {0 -¥ - 1> ; 

Map[xl#&« BiPolarBiteeq] // FUtten] 
CDMAEncodedSeq = CCMAEncode [ {-1, 1, 1, -X). Qoldfla^liflt / / taot] ; 



■ CDMA decoding of single Symbol 

The modular output associated with { 1 ) 

ModOutputPluflOfte= Modulator [tj [ CDMAEncod© r { X> , (Goldseqliet // Last) /. (O -> -l>i „ 
NumberOfCurvea -* 2, ModulatingPulae OptPulfieScaled, SamplingXnterval -» T/4}/ 

This is a primitive decoder built to study the autocorrelation. This will help in decoding 
Take [ModOutputPlusOnd, 10] 

i0 k e ?llZl? °^ 5 121 32 X ' °- 4312 57 + 0. 748432 I. 0.130196 + 0.363909 1. 
-0.103565 + 0.653405 1, - 0 . 510127 + 0 . 69136 1 , -0.748423 +0 431231 1 
-0430757 - 0*74097 I}' ' 0 ' 853326 " 0 ' 183 682 1 ' —0.69115 - 0.510321 1 § 

Pr±mit:iv©CPMJUl©c©iv©r[ModOutput_ . Golds o<z_. Sample , OverSamplino 1 s = 
MOduln[{xl. xiStACe. xSUpdafce, ac4, xSStata, atfi> # J 

xX ^ Partition [Modoutput , ovarsmnpllnsr] // Transpofld // # [ [Sample] ] & s 
X23tftte = GoldSaq / . {O -> -X}; 

xsstate,= Table[(-l)*Mod[i. 2], {i. 0, Irfuisrth [estate] - 1) ] 
xSupdate zt= Module [{}. 

estate - Rotat©Right[x2Sttttft] ; x4 = FoldList [Plus, 0, estate] //Rest 
* 5StatB = RotateRiffiitCxSState]/ xl (xSstate x4> // Apply [Plus, 
Table r*3Updafca, {i, i, Z.engtl»tGoldS©<a]>] J 4 J ' 

Wc make it more efficient 

Priaitiv6C3)t!AReceivBr2 [ModOutput. , GoldSe<a_, Sample.., oversanplin? 1 : = 
Moduie[{xX, x2State, x3updato, x4, x5staee. x6). 
xl = Partition[KodOutput # OverSanpling] // Transpose //#[ [Sample] 1 &; 
x3stato - aoldSeq/. {0 -> -l); i l j j , 

f state = Tatole[(lj -Mod[i # 2]. (1, o. r.ena«i[x2 state] ^1)1; 
x3update in Modulo [Q. J /J ' 

xastate ^ RotataRignt [x2 State] j x4 = FoldLiot[Plua, o, estate] // Reat // i--#fc. 

Table [x3TJ*ia a fco, {i, i, Lexi ff th[Gold6©q]>]] J J * 

Tom = PrimitiveCDMWleceivartModOutputPludOne, ((kjldaaqliat // Laat) , 1 # 4]- 
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Graphics 



Tom = Prinw.tiv©CDMAK©coivftr2[Md>dOutimt36iluoOii©, (Coldsedlist // Last) # 1, 4J , 
iii at Plot [Tom // PlotJoinod -> True, PlotRange -> All] 



50 



-50 



-100 



L 



50 100 150 200 250 



-150 
- Graphics - 

Tom3 = PrimitivoCDMAR6caivar2tMoaOutputPlusOne3. (Golds^lifib // # [ [3] ] &) # 1. 4 
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- Graphics - 

Tom4 = PrimitiveCDMMlweWtModOutputPluflOn^, (Gold QeQ li st // # [ [4 ]] &) , 1. 4] : 
LiatPlot[Toni4 // Ra # S-lotOTQined True, PlotRange -> All] 




-10 
-20 
- Graphics - 

T™,5 = P r i J ni t i VQ c DM AK« lca l VQlrr M 0 d 0utpilt p lus o n e3. (eoXdseallat // If t [*)]*) . 1. 4] .- 
tiatPlot [TonS //Be, Plot Joined -> Trua. PloCRaxicje -> All] 
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Length [Tom] 
255 

Taka [Tom, 20] 

{1.75501 + 0.329995 1. 0. 0343421 - 1 . 5B366 I , 1.39349 + 29.6902 1, 153.017- 0.957068 1 
-0.71B377 + 165.563 I, -43.6424-2.04777 1. 2.21031+0.2566361, 5.392+2.186421. ' 
-3.13017 + 6.54167 1, 3.61193 - 2.5452 1. -1.69586 + 5 _ 24B64 I, 3.35077 - 1.96193 1 
4.9340B - 5. 10993 I, -1.90672 - 5.29083 1, 3 . 34843 + 0 . 9149B3 I, -1.24192 + 4.93512 1. 
-3 .90572 - 2 .50766 I, -6.70085 + 2.60649 1, 0.886488 - 4.00636 1. 10.5B07 - 2 0B322 I} 



L±3fcPlot [Tom // Ra // Tales [#. 20]fi. FlotJoined -> True, PlotRangft -> All] 




- Graphics - 



Li B tPlot(Tom y / Im / / Ta)fift[#. 20]fi, PlotJolned -> True, PlotRanga -> All] 




- Graphics - 



Take [Tom, 10] 



{1 n 7 5?2,-rr.°"??2 9 f« 1, 0-0343421 - 1.58366 I, 1.39349 + 29.6902 !. 153.017 - 0. 957068 I 

:§:«2I?l^ISii? 3 x^3:Sil3 4a .V.I*5°rt 7 , 7X - 2 - 21031f0 ' 256e361 - s -»»^-»«». 
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LifltPlot[Ton // Re // Taka[# ( 20] &, PlotJoiiied 



-> True, FlotRan^a -> All] 




- Graphics - 

UgtPIot[Tom// lift // Takotfl, 30]&, PlotJoined _> T jrue, PlotHange -> Air) 




2 -5 7.5 To i^rs-^is r? '.Y = t s 5 

- Graphics - 

LifltPlob [Torn //Aba //T^[# # 20]a, Plot Joined -> T™, PlotRajjffe _> All] 




We try a less favourable sequence 
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»adule[{xl. x2state, x3updat«, x4. xSState, x6) , 

x2Stafc* = Rofcat ©night rx2 St afco] f x4 = FoldLlut [Plus, o, x2SCate] // Root J / i * ««, - 
Table [x3 update, (i, l, Length [Golds**] > ] J J Jr 

Tom4 = P3r±m±tivGCDMM©ceiver[MoaoutptitPlTxaOiie3 # (CSoldsaqlidt # [ [3] ] &) , i# 4] , 
Take [Tom4« 10] 

{14.7675 - 4.24542 1, -0.301874 + 6.61008 I. -5.92939 + 29 6039T -ISA fii in ^ 

fi'SJI! : J: fiSfe-.VSSft- 1 '-" * 

Liatpiot [Tom4 // Ra // Takef*, 20J*. FlotJolnod -> t™©, (lotRange -> All] 
50 r 





LiatPlot[Tom4 // iro // Tafce[# f 2Q]& # piotCToined -> True, PlotRange -> All] 




Graphics 



HodOutputHinuaOne = 
ModulatgrtX,] f CDHREn CO de[{-l}, (Goldseqlist // taafc) / . {0 -> -111, 
NuznberOtCuirvM - 2, Modulatingfrulee opt Pule ©Scaled, sanpli^Intorval -» T /4] ; 



TomHi 



PrimitiveCDMAReceiverMinuB [ Wo clout pufcMlmifi One, (Goldaegli 



at // Laat) . 1, 4] / 
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HatPlotfToxnHl // im. Plot Joined. -> True, PlotRaage -> All] 




25 

-25 
-50 
-75 
-100 
^125 
-150 



- Graphics - 

ToxnMl = PrindtiveCDMAReceiver [ModOutputMiruisOne, (Goldeeqliat // taat), l, 4], 
LietPlot [Tomial // In, Plot Joined -> True, Plot Rang© -> All] 




- Graphics - 

Length [ Tom ] 
255 

Talc a [Tom, 20] 

{X 'n%\l\^°'^ll 3 ll^% 1 " 58366 1.39349 + 29.6902 I. 153.017 - 0.957066 1. 

"?*Tini^ 7 f I5 3 ^ X *^"5? ' 642 4 - 2.04777 I. 2.21031 + 0.256636 I, 5.392 + 2.18642 1. 

+ 6-54167 I, 3.61193 - 2.5452 1. -1 . 685B6 T 5 .24064 I, 3.35077 - 1 96193 1 

^loS^'^^n^B^- 1 ^^!.: 5 - 2 ? 083 3 .34843 . 0.914983 i. - 1 24192 + 4 93 512 I , 
-3-905.72 - 2. S076B I, -6.70085 + 2.60649 1, 0.8B648B - 4.00636 1 , 10.5807 - 2.08322 1} 

Checking the correlation properties of the Training Sequence in GSM 

GSM = (0, 0, 1. 0, 0. Xs 0, 1, 1. 1. 0. 0, 0, 0, 1, 0, O, o, 1, 0, O, X, O, 1, 1, 1) 

{0. 0, 1, 0. 0, 1, 0. 1, l, 1. 0, 0. 0. 0, 1, 0, 0, 0, 1, 0, 0, 1. 0. 1. 1. 1) 
GSMae* = c-1- -1. l. ^-1, -1, 1, -1, 1, X, X, -1, -1, -1, -1, i, -i, „x, _x, Xs 

£ if -i^'i^irV)" 1, x ' _1, lf 1 * - 1 * - 1, _1, 1< - 1, _i * i - 
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Autocorrelationsequence [QSMaaq] 

{26, -2. -2. 2. -2. -2. -2, 6, -10, 2, 10. -2, 
-2. -2. -2, 2. -2, -7} 



-2. -2, 10. 2. -10, 6. 



ModOutputQSMs Modulator [L] [ GSHsoq, ftunfa»0£curves - 2, 

ModulatingPulaa -» OptPulaaScalad, SajaplzLn^ntorva.! T / 4 ] j 

TomGSMl = PriaitiveCDKAReCfliver [MoaOUtputGSM, GSMeeq, 1, 4]; 
LiatPlot[TomGSMl // Abs, Plofctfolned -> True, PlotRange -> All] 




10 15V^/ 20 v 25 

- Graphics - 

PrimitiveCDMAR6caiveiir2 [ModOutputGSM, GSMeecj, l, 4]; 

ListPloe [%90 //Abs, PlotJOinad -> Tru B , PlotRange -> All] 




- Graphics - 

PrIinifciTreCDMzUlQceivear [ModOutpntGSM, GSMeaq, 1, 4]; 
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14 



LietPXot[%92 // Abe, Plot Joined -> True, PlotRang* -> All] 



14 

r 2 
r° 
A / 8 






10 




20 


25 



- Graphics » 

^J^^*^^ ' GolclSeq_. S^le., Overs axnpl in* J 

Module[{xl, xastato, x3update, xSstate, x5. x6), J 

xi = Parti tio»[ModOutput, OverSampling] // Transpose // # [ [sanplA] ] & ; 
x2etate =. aoldae*/. {0 -> ^ 

xSState = Table[(-1) -Mod[i, 2], [i, 0. Length [x2 state] - 1)1 - 
xBUpdate = ,= Module [{), x2State - Rot at ©Right [x2 St ate] ; 
X4 = PoldLidt [Plus, 0 r x2State] / / Reet; 
x5 = Join[{l}, x2 state // Drop [#, -1]*]* 
scfi = FoldJJisfc [Plus, 0, x5] //Heat / / I^ftt; 
KSStato = Aot a t e Riarhfc[x5State] # - xl ( xSState x6j // Apply [Plus, 

TomGSMgacondP = Pri«itiv-©CDMAR©c©iv©rGSM2Pulo© [ModOutputGSM, GSMse*. 1, 4] ; 
LiritPiot[TomGSMSecon«ss» // Abs, Plot Joined -> True, Plot Range -> All] 




- Graphics - 

Pr±mltiveCDMAReceivBrGSM2PulaeEfficient[ 

ModC>utput_ , QoldSdQ_, Sample^, Overflampling,] : = 
Module [{xl, x2State, x3U*date, x4, x5state, x5. x6}. 

SaUt^^TJiJ < Ti e >?" mPliag] " // # [ [Sample] ] m, 

= T t b 5", C( rl* ,teaCi * 23 '' {± ' °' ^3t:h[x2state] -1)]; 
x3Update := Module [{} , x2State = RotateRight [x2State] ; 
x4 = FoldLiat [Plus, 0, x2State] //Rest; 
xS = Joln[{lJ, x23tata // Drop[#, -l]fc], 
x6 " PoldListfPlus, 0, x5] //Rest// 

TomGSMSecondPE f f 2 = 
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- CDMA Decoding of Several Symbols with Training sequence Receiver 

First we generate the modulator output. For simplicity we will use a very short training sequence . Let the training 
sequence one of GSM training sequences Let the guard sequences be [ 1.1.1} . Let the data symbols be generated by a 
random 

datal a {1, 1, 1, X, X, 1, 0. 1, 0, 3-, O, X, 1, 1, 0, X, X, 0„ 1, 0} ? 
<iata2 = {0, o, X, 1, 1, O, 1, O, X, 1, 1, X, 0, X, 0* 0, X, X}; 
guard = {X, 1. 1} 

(1. 1. i> 

training = {0, 0, X, 0, 0 # X, 0. 1, 1, X, 0, O r O, 0, 1, 0, O „ 0. 1, 0, 0, X, O, X, X, X); 

The GSM Training sequence is {0.0, 1 .0,0,1 ,0, 1,1,1, 0,0.0,0, 1 ,0,0.0,1.0,0,1, 0„1J ,1 )- In fact tiny short m-sequence can be 
used to characterise the output. 

Only at this point that we differentially encode using the gsm scheme 

frame = Join[guar<3, Oatal, training, data2 r guard] 

General: : spelll : 

Possible spelling orror- new symbol name "frame" is similar no existing symbol "Frame" - 

{1, 1. 1, 1, 1. 1. 1, 1. 1. 0, 1, 0, 1, 0, 1, 1. X, D, 1, 1, 0, X, 0. 0, 0, 1, 0, 

0, 1, 0. 1, 1. 1.0, 0, 0, 0,X,0, 0. 0, 1, 0, 0, X. 0. 1, 1. 1, 0. 0, 1. 0, 1, 

1, 0, X. 0, 1, X, 0, 1, 1, 0, X, 0, 0. X, 1. X, X, X) 

GSMDiffEncod«dPrame[frame_] : = ModuXe [ {xl, x2, x3>, xl = Partition [f rune, 2, X]; 
*2 = Map[Mod[#[[i]] + #C(2]]. 2]&, xl] // or Q in[ {frame [ [X] ]> , 
x3 = x2 /. {O -> 1, 1 -> -X)] 

General: ; spelll : 

Possible spoiling error: new symbol nftme "frame- is similar to existing symbol "Frame". 
frameEncodecl = 0$M£)i £ fBncodedFr ane [frame 1 

{-1, X. 1. X. X, 1, 1, X. 1, -1. -X, -X, -1, -X, -1. 1, 1, -1, -X. 1, -1. -1, -1. 1, 1. 
-1, -X. X, -1, -X, -1. X, 1, -1, 1. 1, 1, -X. -X, 1, 1, -1, -1. X. -X. -1, -1, 
1, X, -X, X. -X. -X. -1, 1, -x, -X. -X. -X. X, -1, -1. 1, -X. -X. -1. 1, -1, X, X. X. X] 

ModOutputFramo3 = Modulator £L) [ CDMAEttCode [f rajnefcncadod, QoldaogliBt // Last] , 

NUmbarOfCurveB -+ 2, ModulfttingPulae h- OptPul^eSCAled, Samplinjlntorval T / 4 J ; 
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CirrfreftOOl = Tabl B [E* (I 2 Pi O.OQ1 j) // N, { j , O, Lensjthi [ModOutputPram©3 ] -l}]i 
CarrOOl = MOdoucpuCFrame3 CarrfreqOOl; 
( Save["C^j"001,ia h / CarrOOl]; 

Carrfreq0005 = Table[E» (I 2 Pi 0.0005 j) // Vt M (d. O, L«agth[ModoutputP ra m B 3 ] -1}]; 
GarrOOOS - ModOutputyramoS carr£rdQ0 005; 
Sava["Carr0005.M", carrOOOS] ; 

Carrfx©«002 s Table [e * (I 2 Pi O . 002 j ) / / {J, O, *engtH[WodOufcpufcF2ratae3] -1}]; 

CarrOOa = MadOutputPrajne3 Carr£reqQ02; 

flavo["Carr002.m" f CarrQ02]; 

uodoutputunEiLGodAdFra&M 3 - 
Modulator[I,] [ CDMABncodd[£ran» /. O ->-l, Goldseijliflfc // # [ (3 ] ] , 
Numberofcurvee -» 2, Modulat insrFtil sa — OptPaleaScaled, saicplin&liitQinr^i t / d] ? 

Sav6 E "ModOiitputUnExxcoaedyrameGSMl.ilceS .m*, UodOutputFramo] 

< < MOcLOiitputFi:aiiiBEiiccMlodlGSML.i.lcQ . m; 

Length [ModOutputFrama] 

73440 

AS*cooo5 = Tajc© [Carr0005, {50, 7300}]; 
AFCOO05 // {Ra[#] , Xm[#j ) & // Tranopos© // 

Li8CPlot[#. PlotJoined -> True, Plot Range -> All, AspectRatio -> 1] &j 
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AFC0005 // Taka[# # 200]& // {Ra[ft] . Ia[#])& // TriUflpose / / 
LietPlat [ft, PlotJoined -> Tnio, PlotHaage - > All, Aspect. Ratio -> 1] &; 




In the PrimitiveCDMA receiver we need to specify ihc sample. In the CDMA synchroniser we discover the sample. 

CrjMACoarseSynclxroniBQx-NBW[MCidOiit:put:_ 0 Gold3eg . Throsbold.* OvarSampling^) -- 
Module [ {xl, v2, x3, x4Plue« x4Minus # xSState, 

xficount, x^MaxCorr, ^eq, 3c7Update, x8, x9 ( xio, xll, xl2« xl3), 
xl = ModOutpufc ; 
x2 = aoldSeq / , 0 -> -i; 
3t3 - Codinsrrran3f oznNew[Z>] [GoldBeq] / 
X4P1U9 = x3t[l]] ; 
xdHinufl ' = X3 [ [ 2 ] ] ; 

xSBtate s 7<i]ca[xl, ( Length [x2] «*- 1) Over Sampling] J 
xficount s 1 j 
3c6wavCorr t= Oj 

oocj = Drop [jcI, OverSampling ( X*eng*tli(iE2 ] ♦ 1) ] ; 
x7Update : ss Module [ { ) , 
xB ~ Partition [x53fc ate, OverSaapXing] / / Transpose; 
x9 = Map [ {Drop [ft. -1] . x4Plua, Drop [ft. 1] . x4Minue)& # yfl); 
xlO = Map [ 

{ Afaarim[#£[inn . Aba [Re [ft [ [ 1] ] ] ] . Abfi [Re [# [ [2 ] ] ] ] # Abo [Im [ft [ [2 ] ] ] ] >&, xD] ; 
xll = If [Max[xlO] > Threshold* Throw [ {xlO , ac 6 Count, True}], 

{xfiCmitit:, xlO, xSMavCorr = Max[x£H&xcorr, xlO] , False)]; 
jc 6 Count = x6Count + 1; 

xSState = Join [Drop [x5 State, Over 3 amp ling] 9 TaJce[secx# C*vor Sampling] ] ; 
aeq « Drop [sag;, OvarS&mpliag] S 
xll ] z 

xl2 = Catch[Table [x7Updato, (1. 1* I*ength [ eecj] / Ov-er samp liner > ] ] ; 
Xl3 = X£ [Last [jc12 J =a= True, 

CDMAFinoBynchroni oer [ModOutput , afl2 [ [2] ] # x3 , OvarSampling] , 
{"Failed to Coarse Synchronifle' , False}]) 

Tom 4 o 

CDMACoarsoSynchroni eerNew [ AFC 0005 // Drop [ft, 250 4]&, Goldfloqllse // Laet , 50, 4] 

DeModulatorf 

{{{16.252 - 103. 911 1, 9.19889 6. 9303B 1, -3.61972 + 15.27 63 I, 13 . 0046 - 9 . 41072 1 } . 
{-102.686- 22.7446 1, 7 . 4 9431 - B . 7 45 5 8 1 , 15.0189 + 4. 57178 I, - B . 575 5 9 - 13 . 5 69A X } , 
'{-29. 1474 + 101 . 055 Z. -8.25775 - 8.02866 1, 5.50581 - 14.7022 I, - 14 . 0815 + 7 . 70661 I) , 
{99.0255 + 35.4352 1. -8.53133 + 7.73733 1. -14.3275 — 6.4181 I ( 6.B0721 + 14.5376 I], 
{41.5832 - 96.6051 1. 7.18638 + 9.00033 1, -7.30507 + 13 .8962 1, 14.9364 - 5.Q8096 I) ) , 
{{24.6533 - 96.8492 1, 10.1266 + 6.85283 1, 1.46802 - 2.37691 1, 12.50B7 - 10.488 I), 
{-95.1101 - 30.6858 1, 7.47516 - 9.6763 I, 
-2.28004 - 1.61437 1, -9.68188 - 13 .1425 1}. [-36.5973 + 92.9956 1, 
-9.16784 - 8.06799 1. -1.75435 + 2.17418 1, -13. 7245 + B. 83755 I}, 
{90.5141+42.3643 1. -8.62898+8,66312 1, 2-05973+1.8874 1, 7.95834*14.2523 1}, 
{47.9641-87.6755 1, 8.1042+9.155911, 2.01301-1-937151. 14.7239-7.047721}}}] 
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TomQOOS - 

CDMACoarfia3yllchroniserN©w[CarrQ 0 0 5 // Drop[ft ( 250 GoldflAglist // Last « 100. 4]; 

CDMAFineSyncnroxxi ser [ModOutput:_, 
f TnrecnoldC©rrelaticiCoun.fc_, Coyr-e3.atln.ir3 bp - Qversampling_] : - 

v Module [{xl, x2, x3, x4, x5, acff, Beq # x7tfpdatte, xSFirat, 

xSSocond, xSPirat, at 9 Second, xlOFirsb, xlDSecond, xll, xl2, xl3, xl4) r 
xl = X£ [Tb^rsalioldcorrelatioCouxit > 3, 

ThraohDldCorrelatioCount -3, ThxaeholdCorrelatioCount] ; 
x2 ±± covrelatingSeq; 

x3 = Drop [ModOutput, xl QverSanpling] ; 
ac4 = CDMAPoeitionFinder [jc3 , x2 , OVArB&npling] ; 
x5 CDMAPoaitionrinder [Drop [x3 . Length [x2 ] OverSajnplixi&] » x2, OverSanpllng] ; 

x6 = CDMAFositionFinder [Drop[x3, 2 Length [ac2 ] Over sampling] , x2, OverSampling] ; 
X7 = PofiitionAverager [ {x4 [ [ 1] ] f xS [ [1] ] , x6[[l)]}] ; 
xSFlrsfc = Drop(x3, (x7 [ [1] ] - 1) OverSainpXing] // 

Partition [ft. Over Sampling] & // Transpose J/ #[ [x7 [ [2] ] ] ]&; 
xSSocond = Drop [jc3 , x7[[l]] OverSaaiplipg] 

Partition [#, OvorSanpling] & // Trandpoae // # [ [x7 [ [2] ] ] ]&; 
x9Firet = xBPirat // Partitioa[# / Length [:*2 [ [1] ]]] 
xlO First - Map [Function [x, Map[x. #&, x2]], xSFirst] ; 

x9 Second = x8 Second // Parti tion[#, Length [x2 [ [ 1] ] ] ] & ; 
aclO Second c Hap [Function [5c, Hap[x. x2] ] , xSSeconfl] ; 

LeModulator [ {xlOFirst, xlOSeccmd}] ] 

CDMA^OoitionFindeir [ModOUtput_, Corralatinaflaa - Over£ampling_] : = 

Module [{xSStata, xSCount, a^q, x7 Update, xfl, x9, xlQ. xll, x!2. xl3>, 
xSStete = Take[MoflOutput, ( Length[CorrolAtingE8ij[ [1] ] ] +1) OverSampling] ; 
x 6 Count =1; 

o©« = Drop [ModOutput , Oversampling ( Length [CoriTAlatlsgseg [ [1] ] ] + X) ] ; 
x7update x= Module [{) , 
xB = Partition[x53tata, OverSampling] // Tranflpoee; 

x9 = Map((Drap[# t -1] ♦ CorrelatingSaq( [1] J t Drop[ft, -1] . Comlatrng3 9 q[ [2 J ] , 

Drop[#* 1] . corrolatingSeq[[i]], Drop[#, 1] . CorrelatingSaq[ [2] xBl ; 

x 6 Count - x 6 Count + 1; 

x5State = Join [Drop [x5state, Over Sampling] , Take [oeg , Over Sampling] ] ; 
seq = Drop (aaq. Over Sampling] ; 

xlO = Tabla [^Update, {1, 1 # 10}]; 

anil = Mapindexad [wax [ (Abs [Re [#] ] , Aba [Xm [ft] 1 > 3 &- xlO. {3}]; 
xl2 c Hapindexed [Apply [Plus , #] t, xll, {2}]j 
xl3 ^ Position [xl 2, Max[xl2] ] ] 

Wc need lo define a position averager. We for now just take the first element 



POBiCionAvaragar[PoflitionLi8t_] r= First [Positioning t] 
Tom4 = 

CDMACoaraefiynchroniserNew^TaatData // Drop [ft, 250 4]fc„ Goldpeqliat // Last r loo, 
DeModulator [ 

({{-157.056 -r 0.691379 I. 7.37139 - 19.8305 1. 41.1322 + 1.946931. - 16 . 2B4 3 - 23.4715 I) 
{-0.691379- 157.056 1. 18.8305 + 7.37139 1. -1.94693 + 41.1322 1. 23.4715 - 16.2843 1) 
{157.056 - 0.691379 1. -7 . 37139 + IB . B305 I . -41.1322 - 1.94693 1. 16.2043*23.4715 1) 
{0.691379 + 157.056 1. -IB. 8305 - 7.37139 I. 
1.94693 - 41.1322 X, -23.4715 + 16.2843 1}. {-157.056 + 0.691379 1, 
7.37139 ^ 18.8305 1, 41.1322 + 1.94693 1. - 16 . 2B43 - 23 . 4715 1} } , 
{{-169.734 - 0.510132 1. 6.8B71 - 20. 5026 I, 6 . 2 A 607 + 1 . 479 47 I . -17.4944 - 21.5101 1), 
{0.510132 - 169.734 1, 20.5026 -r 6.B871 I. -1.47947 + 6-24607 1, 21.5101 - 17.4944 1). 
{169.734 + 0.510132 1. -6.9071 + 20.5026 1, -6.24607 - 1.47947 1, 17.4944 + 21 5101 I) 
{-0.510132 + 169.734 1. - 2 0 . 5026 - 6 . 6 971 X . 1.47947 - 6.24607 1. -21.5101 r 17. 4944 I) 
{^169.734 - 0.510132 1, 6.B871 - 20. 5026 I, 6.24607 + 1.47947 1 
-17.4944 - 21.5101 1}))] 



11-03-98 16:41 00 44 1276 677720 

MAR. ' 98 (WED) 16 = 47 -NMP PATENTS UK 

CDMAReceiverFoi9rCo€!eDiffEnCOdeSiudicsAFCPattni.nb 



P 35 R-780 Job-222 

FAX:00 44 1276 677720 



Tom4 [ [1] ] // TranapoBQ 



{ { {-157.056 + 0.691379 I 

{-169.734 - 0.510132 I 
{{-0.691379 - 157. 0S6 I 

{0.510132 - 169.734 I, 
{ {157 .05 6 - 0. 69137 9 I, 

{169.734 + 0.510132 I , 
{ {0.691379 + 157 .056 1 , 

{-0.510132 + 169.734 I 
{{-157.056 ♦ 0.691379 I 

{-169.734 - 0.510132 I 



. 7.37139 - 18.8305 I 
, 6. 9871 - 20.5026 I, 
, 1B.S305 + 7.37139 I 
20.5026 + 6.6871 I, - 
-7.37139 + 18 . 8305 I 
-6.B871 + 20.5026 I, 
-0.0 . 9305 - 7. 37139 I ( 
, -20.5026 - 6.8871 I 
, 7.37139 - 18.0305 I , 
, 6. 0B71 - 20.5026 X, 



, 41.1322 + 
6.24607 + 1 
, -1.94693 
1.47947 + 6 
, -41.1322 - 
-6.24607 - 
. 1.94693 - 
, 1 .47947 - 
41 . 1322 + 
6.24607 + 1 



1.94693 1. -16.2843 -23.4715 1} 
.47947 1, -17.4944 - 21.5101 I}} \ 
+ 41.1322 1. 23.4715 - 16.2043 1} 
24607 1, 21.5101- 17.4944 1}}, 
- 1.94693 1, 16.2843 + 23.4715 1} 
1.479471, 17.4944+21.51011)}, 
41.1322 1, -23 .4715 + 16.2843 1} # 
6.24607 X, -21.5101 + 17,4944 1}} 
1.946931, -16.2843-23.4715 1} 
.479471, -17.4944-21.51011}}} 



Take[Tom5[£lJ ] // Trnnnposa, <8, 10}] 



t( {J-fI i?-r 0 - 929033 -7.5712B + 19.051B I, -40.9109 - 1.74704 1. 16.522 + 23 6653 1} 
{169.613 - 0.714545 I, -7.42576 + 21.3328 I, -5.4158 - 0.94077B i .18 7191 + 21 5891 t \ i 
^ll'llll " Xt -1-02745 + 156.682 X, 23.9098 + 15.9863 X. - 1.57069 ■ L 41 4fi52 II } ' ' 

/f 2 ?"!^?*- 6 :! 6 !^ 1 ' °-47B141^169.715l, 21.5305 + 17.4644 1,-1.46045 - 6 27B06J}} 
{i rl-iln, 9 ^i 8 cS?2 5 Z ',1V 056 + °-69l379I, 16 . 2B43 - 23 . 47 15 I, - 41 . 13 22 + 1 94693 I > 
{-6. 8B71 - 20.5026 I. 169.734 - 0.510132 I, 17.4944 - 21.5101 1, - 6 . 24607 + 1 47 947 if } } 
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Map [Aba, Tom6[[l, 2]]] 



{{99.9069, 
{99.9377, 
{99,9377. 
{99.9377, 
{10.9504. 
(10.9205, 
{10.9205. 
{99.9069, 
{10.9504, 
{101 . 129. 
{10.9205, 
{99.9069. 
{10.9504. 
{101.129, 
{10.9205, 
{99 . 9069. 
{11.1392. 
{99 -9377, 
{10 . 9504. 
{99.9069, 
(10 . 9504, 
{101.129. 
{10.9205, 
{99-9069. 
{11 . 1392, 
f 10 .9504, 
{11.1111. 
{10.9504, 
{10.9205, 
{101. 129. 
{11.1111, 
{10.9504. 
{11 . 1111, 
{11.1392, 
{99.9377. 



12 .1968 
12.2274 
12 .2274 
12.2274 
111. 653 
111-621 
111 .621 
12.1968 
111.653 
12 . 0176 



17.2325. 16.2024} ( 
17 .2666, 16 .1692} , 
17.2666, 16.1692} , 
2.82297, 16.3595}, 
17.2325. 16.2024}, 
2 .74759. 17.0524} , 
111.621, 17.2666. 16.1692}, 
12.196B, 2.B2297, 16.35B5}, 
17.2325. 
2.74759. 
17 .2666. 
2.82297, 
16.5597. 



2.82297, 16.35Q5}. {99.9377, 12.2274. 2.7937. 16.3237}. 
2.7937, 16.3237}. {99.9377. 12.2274, 2.7937. 16.3237}. 
2.7937. 16.3237}, {99.9377, 12.2274. 2.7937, 16.3237}, 
2 JL 7 ???1 16.3237}. {101.16, 12.0496. 2.71211. 17.0171}, 

{10-9205. 111.621. 17.2666, 16.1692}. 
{10.9205. 111.621, 17.2666. 16.1692}, 
{11.1111. 110.39, 16.592, 16.2327}, 
{101.16, 12.0496. 2.71211. 17.0171}. 
{11.1111. 110.39. 16.592, 16.2327}, 
{10.9504, 111.653, 17.2325. 16.2024}. 
{11.1111, 110.39, 16.592. 16.2327}. 
{101.16, 12.0496, 2.71211. 17.0171}, 
16.2024}, {11.1111, 110.39, 16-592, 16.2327}, 
17.0524}, {10.9504, 111.653, 17.2325. 16.2024} 
16.1692}. (11.1111, 110.39, 16.592. 16.2327}/ 
16.3595}, {101.16. 12.0496, 2.71211, 17.0171}. 
16.266B}, {99.9069, 12.1968, 2.82297, 16.3585} 
2.7937. 16.3237}. (101.16, 12.0496, 2.71211, 17.0171}, 
17.2325, 16.2024}, {11.1111, 110.39, 16.592. 16.2327}, 
{101.16. 12.0496. 2.71211, 17.0171}, 
(11.1111. 110.39, 16.592. 16.2327}, 
{10.9504. 111.653, 17.2325. 16.2024}. 
{11.1111. 110.39. 16.592. 16.2327}, 
,„ 2.B2297, 16.3585}, (101.16, 12.0496, 2.71211, 17.0171}, 
H?"™' J- 6 -" 67 ' 1$. 2668}. {101.129. 12.0176. 2.74759. 17.0524}, 
111.653. 17.2325. 16.2024}. (10.9205. 111.621. 17.2666, 16 1692) 

16 - 592 ' 16.2327}, {101.129. 12.0176. 2.74759. 17.0524}. 
111.653, 17.2325, 16.2024}, {10.9205, 111.621, 17.2666. 16.1692} 
111.621, 17.2666, 16.1692}. {11.1111, 110.39. 16.592. 16 2327V 
12.0176, 2.74759, 17.0524}, {10.9504, 111.653, 17.2325, 16.2024} 
110.39. 16.592, 16.2327}, {101.129, 12.0176, 2.74759 17 0524} 
Hi*SI 3 *-.i 7 «2 25 ^ 16 - 2 " a >' < 10 -9205, H1.621. 17.2666. 16.1692}, 

16 A 5 ? 2 ' 16.2327}, {101.129. 12.0176, 2.74759, 17.0524}, 
110.422, 16.S5B7. 16.2668}. {99.9069. 12.1968. 2 B2297, 16 3585} 
12.2274, 2.7937. 16.3237}, {99-9377, 12.2274, 2.7937, 16.3237}} 



111.653, 
12 .0176, 
111.621, 
12 .1968. 
110. 422, 
12 .2274, 
111. 653. 
12.1968, 
111.653, 
12 . 0176, 
111.621, 
12.1960, 



2.82297. 16.3585}, 
17.2325, 16.2024}, 
2 .74759, 17 .0524} , 
17.2666. 16.1692} , 



H&X0005 = 

T^cif ^Y/ 3333 //A **' I -C # t[3]]] //JU> S . »[#[[«]]] //JU> B , I»[#[[4]J]V/ia,B}]>i; 



{{104.351, 15 
(101.055, 14 
{93.6035, 15 
{67.0266, 22 
{86.7967. 16 
{98.9714, 14 
{104.302. 15 
{103 .704, 15 
{96.1447, 15 
{S3.22BB, 15 
{62.9341, 16. 
{96.5637, 15. 
{103.567. 15 
{91.0436. 27. 
{83 .8748, 27. 
(86.6792, 16. 
{7419529, 21. 
{B6 . 0108, 22 . 



1379] , (103 .911, 
7022}, {99.0255, 
(90.7175, 
(70. 8Q01, 
{90 .7238, 
{91. 027B, 
(105 .346. 
{88.1543, 
{93.9076, 
{79.0562, 
{87 . 1763, 
(99-7552, 
2264}, {104.931, 
8733}. {89.5104, 
3795}. {91.0829. 
337B}. {67.5876. 
37} . {77 .7339, 19 
4409}, {96.2752, 



2762} , 
9468} 
2927} , 
9013}, 
2963} , 
2097} , 
3897} , 
9297} , 
4519} , 
3294} 



15.2763}, {102.686. 15.01B9} 
14.5376}. {96.6051, 14.9364} 
15.6415}, (71.4359, 25.1654} 
21.9S11}, {74.4539, 20.868B} 
15.4739}, (86.1537, 22.7114} 
25.3506}, {91.S367, 25.7331) 
15.8922}, {91.9519, 27.7219} 
27.846}, {85.9746, 27.2075}. 
15.3B03}. {74.6965. 25.B353} 

16.0518}, {71.3823, 22.4556}, {74.383^ 20^923} 
15.7147}, {B3.B39, 21.6473}. {B5.5445. 22.25} 
24.5667}, {90.6053, 25.0753}. (91.3309, 25.5023} 
15.891}, (92.2522, 27.2697}, {104.992. 15.4722}. 
27.4649}, {B7.6329, 27.2718}, {100.78, 14.7073} 
26.501}, {78.152. 25.9901}, (74.6068. 25.3375} 
24.357}, {66.9937, 22.9638), (74.5498, 16.5852}, 



{67.5603. 23.B519} , 
{77.6444. 19.7766} , 
{B7.5152, 23.2037}, 
{91.8B71, 26.0709}, 
{90.4855, 27 .2999} , 
{B3 .3983, 26. 8136} , 
(86.612, 16.2707) . 



.7041), {BO. 5993 , 20.0B31}, {90.616, 16.0026}, 
15.4329), {98.8878, 15.1306}, {100.942. 15.3636}} 
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ListPlot[Max0005 // Transpooo // PlotCroined -> Trua] 



( 



100 




LiatPXot [MaxOOOS // Transpose // #[ [2] ]&, PlatJoined ^> True] 

28r 



/l 




v 1 3 ^ 



40 




t£0 



Ma»[{Max[{RaI#[[l]3] // Atoa, //fbs, Re [# [ [2] ] ] // Xba, im[#[[2]]l //*bs}l, 

M^[{Ho[#[[3]]] //Atofi, Ita[#[[33]] //Al,*, R©[#[[4]]] // Atos, Im[#I[4]]] //AfaB>])&, 
TomS [ [1, 2] ] ] 



{{1*9. 
{169 , 
{169 
{169, 
{169 , 
{169 . 
{169, 
(169. 
{169 , 
{169. 
{169. 
{169 . 
{169, 
{169. 
{169. 
{169 . 
{169. 
{169. 



715, 21 
734. 21 
734. 21 
734. 21 
715, 21 
794, 21 
715, 21 
794, 21 
715. 21 
734, 21 
715. 21 
794, 21 
715, 21 
715, 21 
715, 21 
794, 21 
715, 21 
794, 21 



. 5305} , 
.5101} , 
.5101} , 
. 5101} , 
. 5305} , 
-6095}, 
.5305}, 
.6095) ( 
. 5305) , 
. 5101} , 
. 5305) , 
-6095} , 
-5305}, 
. 5305} , 
.5305} . 
.6095}, 
.5305} , 
.6095} . 



{169 

{169. 

{169 

{169. 

{169. 

{169. 

(169. 

{169. 

{169. 

{169. 

{169. 

{169. 

{169 . 

{169 . 

(169. 

{169. 

{169. 

{169. 



734, 21. 
734. 21. 
B13, 21. 
734. 21. 

ai3, 2i. 

715, 21. 
813. 21. 
715, 21. 
313, 21. 
813. 21. 
813 , 21 . 
715, 21. 
813. 21. 
734, 21. 
734, 21 . 
715, 21. 
734. 21. 
715, 21. 



5101} , 
5101) , 
5B91} , 
5101} . 
5991} . 
5305} . 
5A91}, 
5305}. 
5891} , 
5891>, 
5691} , 
5305} , 
5891} , 
5101} . 
5101} , 
5305} , 
5101} , 
5305}. 



{169 

{169 

{169 

{169 

{169, 

{169. 

{169 

{169 

{169 

(169. 

{169, 

(169. 

(169 . 

{169. 

{169 . 

{169. 

{169. 

{169. 



734. 21 
734. 21 
715. 21 
734, 21 
715. 21 
734, 21 
715, 21 
734, 21 
794, 21 
715. 21 
T15. 21 
734, 21 
794, 21 
813, 21 
734. 21 
B13, 21 
813. 21 
734, 21 



.5101}, 
.5101} . 
.5305} , 
.5101} , 
.5305} , 
.5101}, 
.5305}. 
.5101) . 
.6095}, 
.5305} . 
. 530.5} , 
.5101} , 
. 6095} , 
. 5891} , 
- 5101} , 
.5091} . 
.5891), 
.5101) . 



{169 
(169 
{169 
(169 
{169 
(169 
{169 
{169 
{169 
{169 
{169 
{169 
{169 
{169 
{169 
{169 



.734, 21 
.813, 21, 
.813. 21. 
.313, 21. 
.813, 21. 
.813, 21. 
-715, 21, 
.B13. 21. 
, B13. 21. 
-G13, 21, 
.794, 21. 
.794, 21. 
-813. 21. 
.794, 21. 
-794. 21. 
.734, 21. 



5101} , 
5B91}, 
5891} , 
5891} , 
5891) , 
5891} , 
5305} . 
5891) , 
5891) , 
5891} . 
6095} . 
6095}. 
5891) , 
6095} , 
6095}, 
5101}} 



damodframe = sign [Re [ TomS [ [1, 1] J Se<jX] } /. {-1 -> 1. x -> 0} 



{1, 1, 1, 1, 1, l. l. l, o, 
1, 0, 1, 1. 1, 0, 0, 0. 0, 
0, 1. 0, 1, 1. o, 1, 1, 0. 



1. 0, 1, 0, 1, 1, 1. i # i # o, 1, 0, O. 0, 
1, 0, 0. 0, 1. 0, 0, 1. 0. 1, 1, 1, 0, 0. 1. 
1. 0, 0. 1, 1. 1. 1) 



0. 0, 

1, 1, 



The first and lasi bits have been lost in the processing 
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22 



r 



£rame 

{1, 1. 1. 1, 1, 1, 1, 1, 1, 0, 1, 0 , 1, 0, 

0, 1, 0, 1. 1. 1. 0, 0, 0. 0. 1. 0. 0, 0, 

1, 0. 1, 0, 1. 1. 0, 1 , 1. 0, 1. 0, 0, 1, 

72 



1, 1, 1, 0, 1, l f 0, 1. 0, 0, 0, 1, 0, 
1. 0, 0, 1, 0. 1. 1, 1, 0, 0, 1, 0, 1, 

1. 1, 1, 1) 



truneatodfrane = frame // Drop[#, i]fc // Drop[#, -1]& 



[1, 1, 1, 1, 1, 
1, 0 f 1, 1, 1, 
0, 1, 0, 1, 1, 



1. 1, 1, 0, 1, 0. 1, 0. 1. 1, 
0, 0, 0, 0, 1, 0. 0, 0. 1, 0, 
0, 1, 1, O, l f 0, 0, 1, 1, 1, 



denodframe 



t rune at ed franc 



1, 0. 1. 1, 0, 
0, 1. 0, 1, 1, 

1} 



1. 0, 
1. o. 



0, 0. 

1, 1, 



{0, 0, 0, 0. 0, 0. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, o, 0, 0. 0. 0, 0 0 o 
0, 0, 0. 0, 0, 0, 0, 0, 0, 0.0. 0. 0, 0. 0, 0, 0, 0,0, o, o. o! o! o! o' o' o' 
0, 0, 0 , 0, 0, 0, 0, 0, 0, O. 0, 0, 0, 0, 0, O) ' * ' 

Tom6 = CM»CoardBSyiicli^aiiiBer[MociOutpiat.OziEiicoaadFxame3 //Drop[*. 250 dU 4 
aoldflaaliiifc // #[[3] Jfc , 100, 4] J 

S Aborted 



Tam6 [ [1, 1] ) Seal > / Re 

{ l" 7 - 37 ' -157.056. -157.056, -157.056. -157.056. -157.056, -157.056. -156 862, 
5.75774, -9.3815. 9.3B15, -9.3B15, 9.3815. -9.18122. -157.37. -156. B62. 9 75774 
T?-; 1 ?" 2 ' o 1 !?^ 76 ' 975774 ' -9-3B15. 9.18122. 157.37. 156.862, - 9 . 7 5774 ," 9 . 1B122 , 
Hl'lZi' S' 3 ? 1 ?! -9-18122. ^157.37. -156.862, 9.55746. 157 . 37 , 157 . 0S6 , 

l 5 f«!P J -9-75774, 9.16122. 157.37, 156.862. -9.75774. 9.16122. 157.176 -9 75774 
9.3015, -9.10122, -157.37, -156.662, 9.55746, 157.176 -9.75774 9 3615 -9 18122 
-157.176. 9.75774, -9.3915, 9.3815, -9.1B122 -157.176. 9!75774 -9 1122 ' 
-157.176. 9.75774, ^9.3815. 9.18122, 157.176. -9.55746, -157.37! -157.056; -157.056} 

Toxn6 [ IX, X] ] Se^I // im 

{ "^2'?™' "I 0 - 38 ". -10.3815, -10.3815. -10.3815, -10.3815. -10.3815, -10.62, 
T^ 5 *™? 9 ' 15 - 6783 ' "15.6783. 15.6783, -15.6783. 15.9004, -10.0454. -10.62 -15 3789 
JS'f 222- 7i°^ 2839 ' "15-3789. 15.6783, -15.9004. 10.0454. 10.62. 15 . 3789 -15 .900 4 
^2 39 ,'c 1 ?^2Z e9 '^ 1 S- 6783 ' l 5 - 9 °04. -10.0454, -10.62, -15.601, 10.0454, 10.3B15, 
10.62, 15.3789. -15-9004, 10.0454, 10.62, 15.3709, -15-9004, 10 2639 
15.37B9, -15.6783, 15.9004, -10.0454, -10-62, -15.601, 10.2939 15 37B9 
-15.67B3, 15.9004, -10.2B39. -15.3789, 15-67Q3. -15-67B3, 15 9004 -10 2839 
-15.3789. 15.9004, -10.2839, -15.3789] 15.67B3! -15.9004 10.2839; 15 601 
-10.0454, -10.3815, -10.3815} 

T©m6£[l, 2]]SACtX // Rd 

^n^P 1 '^ 1 ^,: 73 ,^ - 16 9- 7 ". -169.733. -169.733. -169.733, -169.733. -169.638, 
2 IA'*7 0 ;^ 5 ^A°-^ 5 ;7 10 ' 425 ' 10 - 425 ' -9.B6944. -169.751, -169 . 638, 10 . 444, 
TfiQ Hi ' h^'^^'^i*^ 23 ' 9 - B 6944. 169.751, 169. 63B, -10.444, 9.86944, 

-l Q - 444 « 10-425. -9-06944. -169.751, -169.630, 9.86847, 169.751 169 733 
,t 3 kU B ' ; 1 2^f 44 ' 9 - B6944 - 169-751, 169.630.-10.444, 9 . B 6944 16 9 6567 -10 444 * 
1 ?ca 2 k/ 9 ;5 6 ? 44 ' " 16 9-751, -169.638, 9.B8B47, 169.656, -10-444. 10.425, -9 86944 
~}ll m lH' 10 * 444 ' -10.425, 10.425, -9. 86944, -169-656. 10.444. -9.86944 ** Ub * 44 ' 
-169.656. 10.444, -10.425. 9.86944, 169.656, -9-60B47, -169.751, -169.733 -169 733} 
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Tont6 - CbMACoaraeSyncliroxxlaexr[MoilOut:put:nzi«z3cadailFzraiaa3 //Drot»r#. 250 4 1 
GoldBeqlist: // #[I3]]& , 100. 4] J 

DGModulator[{ {10.3615 - 157.056 I. 157.056 + 10.3B15 I , -10.3815 ^ 157.056 I 

"J". °^6 " 10 ' 3815 X ' 10 -3915 - 157.056 X, 157.056 + 10.3815 I, "-10.3815 + 157. 056 I 
-156. B62 - 10.62 I, 15.601 + 9.55746 1. -157.176 — 10.2039 I, 15.601 + 9 55746 1 
7~ 7 »J 76 ' 10 - 2839 I. 15.601 + 9.55746 1, -157.37 - 10.0454 1, 10.3015 - 157.056 1 
156.B62rl0.62l, - 15 . 601 - 9 . 55746 X . 157.37 + 10.0454 1, -10 62 + 156 862 1 
9 ; 5 " 46 " 15,601 X ' 2839 + 157. 176 1, 9 - 7577 4 - 15 . 37 99 I, -15.6783 - 9.3 Bl 5 I 

-9.18122+15.9004 1, 10.2839-157.176 1, -9.75774+15.3789 1. 15.9004+9 18122 X 
Til 7 i« 6 ;i 0 *5 839 1 - 15. 601 + 9. 5S746 X, -157.37 - 10.0454 1, 10 . 3815 - 1S7 \ 056 I . 
i 5 fi?S2*3-2-S5J' -15. 37H9 - 9. 75774 1. -9.3815 + 15.6783 1. 15.67B3 + 9.3815 I, 
^o^c 5; l h 9 25 41 ' -10. 2B39 + 157.176 1, 9. 75774 - 15.37B9 I, - 15 . 6783 - 9 . 3 B 15 I , 

J?2 + J5"22?2 r ' 10 '?839 - 157.176 1, -9.75774 + 15.37B9 I, 15 . 9004 + 9 . 10 122 i , 
7cf «* 7 6 ~** C zVJ* X. 15.601 * 9.55746 X. -157.37 - 10.0454 1, 10.3B15 - 157.056 I, 

I JL% 2 aZ - 1 S-37B9 - 9.75774 X. -9.18122 + 15.9004 1. 10.2839 - 157. 176 I , 

~^ln fni 6 ii 5 ^?i Z ' 10.0454 - 157.37 1. 156.862 + 10.62 1, -15.601 - 9.55746 1, 
o «i2« + i2"22? 91 ' Ti 5 -; 6 ? 1 " 9 ' 55746 1 ' 15*?. 37 + 10. 0454X. -10.62 + 156.862 1, 

iln " 1 ?« fi Si-?Li -1 0.0454 + 157.37 I, -156.862 - 10.62 X, 15.601 + 9.55746 I, 
~,ll'nl%~**'*lHl' IS. 3709 + 9. 75774 1, 9.18122 - 15.9004 1, -10.0454 + 157.37 1, 
7 if 7 4?f 6 " ^-P- 30151 ' 10.3815 - 157.056 1, 157.056 + 10.3815 1). {-7.65342 - 169 . 73 3 I , 
"!'^3 - 7.65342 I. 7.65342 + 169.733 1, -169.733 ^ 7.65342 1. -7.65342 - 169.733 1 
^l^ll* 7 '. 65342 *' 7.65342 + 169.733 1, - 169 . 638 „ 6 . 41574 X . 17.2252 + 9.06047 1, 
+ 6-44773 I. 17.2252 + 9. BB847 X, -169.656 + 6. 44773 I, 17.2252 + 9.908471 
T^^e? 1 t 7 ; 68541 -7.65342 - 169.733 1, 169.638 - 6.41574 1, - 17 . 2 252 - 9 . 88 847 X 
JS 9 iZ5 ,Z-S25 41 Z ' 5- 41 574 + 169.63B I. 9.88B47 - 17. 2252 I, 6.44773 + 169.656 1, 

?a 117 22^' 10 - 425 x - -9. 86944 +■ 17.25S3 I, -6. 44773 - 169.656 I, 

~J2a 4 25 ^ 16-3B2 I, 17.2553 +9.66944 1, -169.656 + 6.44773 1, 17.2252 + 9.88847 1 
-169.751 + 7.68541 1, -7. 65342 — 169.733 I, 169.638- 6.41574 1, -16 . 362 - 10 444 i 
7« 0 ;?? 5 : 16 * 412 X ' 16-412 + 10.425 1, 9-B6944 - 17.2553 X. 6 . 44773 + 169 . 65 6 Z , 
10.444 - 16.382 1. -16. 412 - 10, 425 I, -9. 86944 + 17.2553 X, -6.44773 - 169 656 1 
~?2a 4 2S.. + 1 5'2 Q 5 T ' 17-2553 + 9. B6944 X, -169.656 + 6.44773 1, 17 . 2 2 52 + 9 . 8 8847 X 
"q lkl\\*Xi*lli\l' "2*!5 342 _1 52-Z 33 X - 169. 638 - 6. 41574 1, - 16 . 38 2 - 10 . 444 i , 

f 44 * M^ 2553 -6.44773 - 169.656 1. -9.68047 + 17.2252 1, 
-7.68541 - 169.751 1, 169.638 - 6.41574 1, - 17 . 22 5 2 - 9 . 8 9 847 X , 

U^IIH' r 1 J"2252 - 9.88B47 I. 169. 751 - 7.6B541 I, 6.41574 + 169.630 1, 
9 ;?5 8 il" 1 7.2252 I. 7.68541 * 169.751 1. -169.638 + 6.41574 1. 17 . 22 52 + 9 . 88847 I 
~?f! "6 + 6.44773 I, 16.382 + 10.444 1, 9.86944 - 17.2553 1, 7.68541 + 169.751 I 
-169.733 + 7.65342 1, -7.65342 - 169.733 1, 169.733 - 7.65342 1}}] 

Tomfi [ [1„ 3,] ] S©« // He 

157.056. 157.056. 157.056. 157.056. 157.056. 156.862, -9.55746, 
T« 7 n^ 6 ' Mf^ 4 !' 15 I- 176 > -9-55746. -157.37, -157.056. -156.862, 9 . 55746 , 157 37 , 
~l'lll ASm - 157 - 176 ' 9-75774, -9.3815. 9.18122. 157.176, - 9 . 7577 4 , 9 - 19122 . 
\lh ^nk'Vnlinf' ~l 5 1' 31 ' - 157 .056. -156. B62, 9.75774, -9.3815, 9 . 3 8 15 . - 9 . 1B122 . 
"^7 ' ^kn*nZV ~^l B il^ 9 k lB ll 2 ' 157 -176, -9.75774, 9.18122, 157.176, ^9.55746, 
:J"- 37 . -157 056, -156.062. 9.75774. -9.16122. -157.176, 9.55746. 157.37. 
156 862. -9.55746. -157.176, 9.55746, 157.37, 156.862, -9.55746, -157 37 
-156.862. 9.55746. 157.176. -9.75774. 9.18122, 157.37. 157.056. 157.056, 157.056} 

Tome [ [X # X] j s©qi // im 

:3°^^ 15 ' -10.3815. -10.3815, -10.3815, -10.3815. -10.62. 

^ 5 o^ 9 ' ^S™ 3 ' - 15 -67B3, 15.6783. -15.6703, 15-9004, -10.0454, -10.62, -15 3789 
Jn*?2?S' :J°^S 39 ' " 15 - 3 7B9, 15.67B3, -15.9004, 10.0454, 10.62, 15 . 37 0 9 / -15 79004 
10.2B39. 15.3789, -15.6783, 15.9004, -10.0454. -10.62, -15.601 10 0454 10 3815 
J2-I2i- 15 - 3789 ' - 15 - 90 °4, 10.0454, 10.62, 15.3789, - 15 . 9004, 10 . 2 03 9 10*3815, 
15.3789, -15.6783. 15.9004, -10.0454, -10.62, -15.601, 10.2939 15 3709 

^■nRS!' - 10 -2839, -15.37B9, 15.6793, -15.6783, IS _ 9004 , ' - 10 . 2 939 , 
~J"S'nZS2' 1 ^ 9 ?2f^ - 10 " 2839 ' -15.3789, 15.6783, -15.9004 10.2839 15.601 
-10.0454. -10.3815. -10.3815} 

Tamfi [ [l, a ] ] seql //Re 

^O 6 !^ 51 *^^^ 7 ^ ;^ 69 * 733 ' - 165 - 733 - -169.733, -169.733. -169.733, -169.638, 
10.444, -10.425, 10.425, -10.425, 10.425. -9.86944, -169.751 -169 638 10 444 
;f o e f! 44 ' - 169 -656. 10.444. -10.425, 9.96944, 169 . 751, 169 630 , -l6 444 9 86944 

-JS* 444 ' 10 ' 42 5. -9.86944. -169.751 -169.638, 9 88947. "^"Si. i" 733 
in itl ' n^kili' 9 ;I! 9 SS- 169 -751. 169.638, -10.444, 9.B6944. 169 . 6 5 6 f - 10 444 ' 
let «fi r T S 221* "^ 9 ;o! :L '^ :L69 ( : 63 ^ 9 - BB B47. 169.656, -10.444, 10.425, -9.86944, 
'^ll'ili' 10.444, -10.425, 10.425, .9.06944, -169.656, 10.444, -9.86944 
-169.656, 10.444. -10.425, 9.96944, 169.656. -9.88847; -169 751. ^lel 733, -169.733} 
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T&blo [CDMAPoeit ionFindor [Btodoutpu tUnEncodBdPrame3 // propfft, 250 4*1 20lfc 
Gdlda e qli B t//#[[3]] ft| 4 ] //mttattZ/Abfl/ZMax, J ' 

{1, 231. 350)] 

( { Vi'H%V ?D^oa 2 'ia 2 ;^ 92 ' 1^-4452. 17.6558, 17.6558, 17.655B. 17.6558, 17.4676, 

\V\V> 'ifi 1B.0BB. 18.088, 19.6342, 19.6342, 19.6342, 19.6342. 19.6342. 

13.457. 16.0791, 16.0791, 16.0791. 16.0791. 19.95, 19.95. 19 95 19 95 19 95 

2ri^ 2 ' 2 i 6 60? 93 l' ft ^^f 4 ^ 2 ^°I 31 ' 25 '" 7 ^ 25.3773. 25'.3773. isVaVV^aV.fio'l. 
?af™t ^o 6 ?^. 10 :! 221 ' 20.5255, 26.9209. 26.9209, 26.9209, 26,9209. 18.9705 
IB. 9705. 18.9705, 20.3921, 22.6392. 22.6392, 22.6392, 22.6392 17 4452 17 65SB 

}Z*£f S 2' 17 - 655B ' 17.4676, 17.4676, 1B.0B9, 18.088, 18.088. IB. 088 
\l cZ%V i!'S?S?' Jl-g" 2 ,'- 19 - 6342 ' 19 6342 ' 13.437; 16-0791. 16.0791. 
H'ZZH' If'SI^' 19.95. 19.95. 19.95. 19.95. 19_95, 18.3042, 16.6993, 
^ '25*?' 23 .0831, 25,3773. 25.3773. 25.3773. 25.3773, 24.601, 24.601 24 601 
IB. 7221, 20.5255, 26.9209, 26.9209, 26.9209. 26.9209, 18.9705, 18.9705. IB 9705} 

MOSC [*] 

26.9209 



Map [Max. *s] 



{17.6558. 
19.6342, 
16.0791, 
16.7542, 
18 .7221, 
18.9705, 



JZ"«5S- 5-2 " £552- "" 5 *' l 7 - 4676 ' 17-4676, 18.088, 18.QBB, IB. 088, 18.088. 
J?-52 42 ' I 9 - 6342 ' 19-6342. 19.6342. 13.457, 16.0791. 16.0791. 

I 9 " 95 ' 19.95, 19.95. 19.95. 19.95, IB. 3042. 16.6993. 
ll'llll' 25 - 3773 ' 25.3773, 25.3773, 25.3773. 24.601, 24 . 601 . 24 . 601 . 
20.5255, 26.9209. 26.9209. 26.9209, 26.9209 IB. 9705. 18 9705 
20.3921, 22.6392. 22.6392. 22.6392. 22.6392, 17.4452) 



References 
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BeglnFackage [ p RiagrunotionflDiffEncodad* » ] 

hinglDivieioa: lueage = "ILlngDivieiori[ModuloRiiig) [ {DAnonlnatoi, Nuaerator) ] o 

{ Quot itant , Ramanidar . Donomiiiator. Fl ag} a. g If xium ^ {1, 0, 0, O, 0,0, 0, O,0 # 1, 
0,0, 0,0. 0,1}. dan => (1,0,0,1,0,1). thon nin9DiviBion[2j [(aen fW m}] = {{l, 
0,0,l,0 # l,l,0,0,a,i),{l,0,i,o,0},{l,o # o,i,o,l),Trua). Whan tha 

tt??Ti tW encounterfl a daviaor the flag la set to False, e.g RingDivieion r 

4] [{{2,1,1,1}, {1,0,0,1}}] {{}. {1,0,0,1}, {2,1, i,i),F A laa) « L 

DropLoadxngZexroo : : usage = 
"DropLeadingZeros [DenominatorGeq] drops tha leading zeros in the sequence, 
e.g DropLeadingZeroa [(3,3,1)] •« (3,3,1). DropLeadingZeroe [{0,0,0}] = {}« 



PooBibleDiviaora 



usage - " PossihleDi visors [HoduloRing] [Order] gives the list 




PolynomlalMuitiplicatiDn: :u^ge - M FolysomlalHuleiplication[HoduloRing1 [ 
polyl,poly2] performs polynomial multiplication moduioAing e.g. 
PolynamialKultiplication[4] [{2,1}, {2,1}) = (!}" 

ModuloMuitplicatiou; : us age = 

-ModuloMUltplication[MofluloRlng] [PrimitivePplynomial] [PQlvl,POlv21 a 
g Mo^lo»Xiltpli ca tionC4][ ( l,2.1,3}Ht3,3,0)!7^ro>] . (il^- 3 

J), {1, 1,2, 2) , (1,2. 3.2,1), (3,3, 1,1,2}, (1,0, 3,1. 31 . f 3-2 1 1 

i'??'??' 1 ' 3 * 5 -*' 3, 3, 0,2), (3,1,3,3,3). (3,2, 2, 1, 1) , {2, 0^ 1^ 

3, 2}, {1,1,2, 2,0), {2, 3, 3, 3,1). (3, x,o, 0,3} ,{3.2,1, X), {3, 1,0, 2,3), {!})- 

CyclicMultiplatlvaGroupz ruaage = "CycUcMultiplativeGroupfModuloRingl f 

Pr^mitxvePolynottial] a.g. CyclicMulti P l a tiveGroup[4 ] [ {1. 2, 1. 3} ] = , / r x \ n 

2l?(^Si > 3 C 2i^ , ' { r^' l 2> .; { ^ 3 ' 3) ' {1 ' 2 ' 1 >> *»* ^ic^itipUtiva^ro;^' 
/ n ;V'?' 2,3}1 = H1>,{1,0,0), (1,0,0,0,0). (1,2,1,0), {2.1.1,2,1), 1 

fl* A A (3.3.1,1,2), {1, 0,3,1, 3}, {3,2,1,1,0}, (1,0.0, 3. 2), 

^'S^^' 0 ' 0 ' 0 ^ {X ' 2 ' X > * (3-2,1.0.0), (1.1,0,1.2). {2, 1,3,1} # 

0), {2,3,3,3.1). {3.1.0,0,3}, {3,2,1,1), {2, 1,0, 2,3}} - 

Zero Sequences: aus&ga = 
-ZeroSwancafltModxiloRi^g] [Order] generate tha zero ida*l ZaroSequencea [41 f 
3]= <<0^.0>*{0*0,2},{0,2,0}.{0.2, Z },{2,0,0},{2,0,2),{ 2# 2,0},{272r2})- 3 1 
z a r 0 p a d:m Bage = -zeroPad[ordar] [Elament] eg Zero*ad[3) [{!)] = {0,0. 1} « 
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ivoAdicExpanoloB! :uaage = ■TwoAdicExpAnfliGn[MoauloRing] [PrimitivePolynomial] e.g- 
TwrciAd±cE3epaiio±oa[4] [{1,2,1,3)] = {{ {O. O, O), {O, O, O) ,£)}.£ (0. O, O) , {O. 0, 1} , 
(2)), {{0,0.0). (0, 1,0), {2,0}}, ££0, D,0), £1,0,0), {2,0, O) ), {{0,0. 0). {2.3. 

1) , {2,2}}, {{O.0,O>. {3. 3. 2). {2, 2,0}), {{0,0,0). {1.3,3}, {2,2,2)), {{0,0,0), 
{1,2.1), {2,0,2)), {{O, 0,1), {0,0.0}, {!)}, {{0,O,l), {0.0.1), {3}), {{0,0,1} , 
(0,l,O), {2,1}). {{0,0,1}, {1,0,0), {2,0,1)). {{0,0,1}, {2,3,1), (2.3)). {{0.0. 
1}* O* 3,2), (2,2,1)). {{0.0.1), {1,3,3), {2,2,3}}. {{0,0,1), {1,2,1}, {2,0, 
3}), {{0.1.0), {0,0,0}, {1,0)}, ££0.1,0), {0.0,1), {1,2}). {{0.1,0}, {0,1, 0), 
{3. 0}),{ £0.1.0), {1,0,0). {2,1,0}), {{0.1.0). {2.3,1), {3,2}). {{0,1.0}, {3.3. 
2}. {2,3,0}}, {{0,1,0} , {1,3.3). {2.3,2)}. {{0,1,0}, {1,2.1}, {2,1,2}), {{1,0, 
O), {0.0.0}. {1,0,0}}, {{l,O,0), {0,0,1}. {1,0,2}). £ (1,0,0), £0.1,0}, {1,2, 

0) }, {{1,0,0}, (1,0.0). £3.0. 0)}, {£1,0,0), {2, 3.1). £1,2, 2}}, £ £1,0.0) ,{3,3, 

2) , £3.2, 0}), £ (1,0.0), £1,3,3), £3.2.2}}, £ £1, 0 , 0 ) , { 1, 2. 1} . {3 . 0 . 2} } , { {2 , 3 . 
1># (0,0,0) , {2,3,1}}. ££2.3.1}. £0.0.1), {2,3,3}} , ££2,3,1), £0,1,0), £2,1, 

1) }. {£2,3,1), (1,0,0), (3,1)), {£2. 3,1}, £2, 3, 3-), (2,1,3) >. £ £2. 3, 1) , £3. 3. 2) . 
{1*1) ># {{2,3,1). {1,3,3). £1,3)), ££2,3,1), (1.2,1), £3,3}), ££3,3,2), (0,u, 
O), £3,3,2)), (£3,3,2) , £0,0,1), £3,3, 0)), ££3,3,2), £0,1,0), £3,1,2}), ££3,3. 

2) , £1,0,0). £1,3,2)). ££3.3,2), {2.3, 1>, £3. 1,0)), ££3, 3, 2), (3,3,2), {1.1, 

2) ), {{3.3,2}, (1,3,3), (1.1.0)). {(3,3.2), (1,2,1), (1.3,0)}, ££1,3, 3), £0,0, 
0} , £1,3*3}), £{1,3,3}, (0,0,1), {1,3.1} ), (£1,3,3), (0,1,0), {1. 1,3)), {{1,3. 

3) , £1,0,0), £3.3.3)), ££1,3,3). £2,3,1), (1,1.1)}, ££1,3.3), £3,3,2), £3,1, 
3})* {{1, 3,3), {1.3, 3), £3.1.1)), ££1.3,3), (1,2,1). £3.3.1}), £{1,2,1}, (0,O, 
0). £1,2,1}), ££1,2,1), (0,0,1), {1,2. 3)}, ££1,2,1), £0,1,0), £1. 0,1)), ££1, 

2.1) . £1.0,0), £3,2,1}), ££1,2,1), £2,3.1). {1.0,3}}, (£1,2,1), £3, 

3.2) , £3,0,1)), ££1,2.1), £1,3,3), £3*0,3)), {(1,2,1). £1,2,1), £3,2,3}} }" 

a: : usage = h cr [MocJuloRing] [PrimitivftPOlynomial] [Twoaaic Element] e.g - 

AutomorphiBraSigma: :u£Jajya = "AutOrmorplii.smsiQrma [Dlod\xlo»in&] [PrimitivoPolynoraial] e.g 
AutcwnorpliloinSiHinar*] [{1.2.1.3}] = £ £ £0, 0, 0) , £ 0, 0 , O] ) , £ £0 . 0, 2 } , £0, 0 , 2) ) . 
££0,2.0}, £2.0,0}}, ££2,0,0), £2,2,0)}, ££0,2.2). £2, 0 , 2 } } , { {2 , 2 , 0 ) , {0,2, 

0) }, £(2,2,2). {0.2.2}}. {{2.0,2}, £2,2,2)), ({0,0, 1), (0.0,1)), {{0,0,3), £0, 
?i 3 > } :JJ°; 2 * 1} ' {2 ' 0 ' :L >>-^ 2 ' 0 - 1 >- C2, 2.1}}, ££0.2, 3), £2, 0.3)}. (£2. 2.1). 
^?'?' :L J } l^^ a '^ 3, • <°' 3 ' 3 >>'H2,0,3},£2,2.3)}, ££0,l,0},£l,0,O)),££O,l, 

2) , £1.0.2}}. ££0,3.0), £3, 0,0}}, ££2,1, O), (3,2, O)}. ££0,3,3), £3,0,2}), ££2. 

^?"?;? } ; { ^ 3 ;2 , J'i^' 2 ' o> ' <:L ' 3 ' 2>) - {{3 * 0 - 0 >' ^-3- 2 >}*<ci^2,2),£i,3, 

2 > h^?i5'2 > 4P^'5 } L' a3 * a ' 2> ' C3 ' !L ' 0> >' {{3 - a ' 2 >- U'1'0)). ££2,3,1), £1, 

?; 1 i > *fP;^ 3> ' <1 ' 2 * 3>) ' {{2 ' 1 ' 1} 't 3 ' 2 ' :L >>'H 0 - 3 * 1 >*t 3 - 0 ' l >>'a2,i,3), 

V'J?'*'° }} * << 3 ' 3 ' D >' (0,1,2}), {{3, 1.2), £2, 1,0)}. {£1,3, 2). (2, 
,',Vl'kP', ±m ° } ' t 2 ' 1 ' 2 )*' £{1,1,2), (0,3,0)), (£1.1.0), £0.3,2}), 

^ 1 ;«'2 > ^ {2 ' 3 ; 2) l'iJ 1 ' 3 ' 3> '^ a ' 3 ' :L >>' { < :L ' 3 ' :L >'< 2 ' 3 - 3 >>'Hi-i* 
? ) ^^?^^ } ^^P;^ 3) ;^^ : !:; :L}3 ' ^^1 ' :L ' 1 >' t0 ' 3 ' 3 J>'^ ^3 - :L ' 3 >-^ 2 • 

r^ } J'iP'^^ ) ' <2 ' :L ' 3)> ' H3 ' 3 ' :L} - 4:0 ' :L - 3 >>'^ :L ' 2 ' 1 >'< :L ' 3 ' 3 >)' 

^^/^? > ^^^^/i } ^K cl ;?'^ } ;P' ^ ' ^}> ' {^3 ' 2 ' l> ' ^3 ' l - 3} >'" l • o * 

3) , £3, 3,1)), {{3, 0.1), £1.1,3)), {£3, 0,3}, £1,1,1} ),£ £3, 2, 3). £3.1.1} )} <• 

unitsRlng: -'uaage = -UnitsRixig[MQduaoRing] [PrimitivePolynomial] e.g. UnitsRincr T4 1 T 
ir'^U^l 1 = {<°*0*^}. (0,0.3). £0,1.0}, (0,1,1), (0,1,2), (0,1,3), £0,2, 1 

1) ,£0,2,3) £0,3,0), £0 3,1} £0,3,2), £0,3.3). £1.0.0}, £1,0,1), £1,0,2}, £1, 

o ^^^^ l ^ ^^i^ 2 ^t l ' :L ' 3 >'< :L ' 2 ' o >'< l ' 2 ' l >'t^2.2).£l.2.3). 

iv'r; J r ^'^ 1} ' (1 ' 3 - 2> * (1 - 3 - 3} * {2 '°' 1} ' ( 2 '0- 3 )* (2,1,0), (2,1.1), £2,1, 

2> AP;^« > ;P;^i > • ^2 ' 2 ' ^> '< 2 ' 3 ' 0 >' t2 ' 3 ' :L >^ 2 ' 3 * 2 >-^ 2 ' 3 - 3 >-t 3 * 

S; h { 2' 0 ? 1} ,' < 3 ' 0 * 2 >'< 3 '°' 3 >'{ 3 '1*0},{3,1,1),(3,1,2),{3.1. 

3), (3,2,0), (3,2.1), £3,2,2), £3,2,3), £3.3,0}, £3,3,1), £3,3,2), £3,3,3) )- 

-Tz :U9^ge = "r[Order_] [x_] b . 9 <r [4 ] [«3 J [ £3 . 0, 3)] = " 
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TracenepraBentation: s us age = "T*ace*iopr*Hexit:atii>n[woauloRiiiff_} [ 

Primitivopolynoiaial_] o-cr- TraceRepr«B©ntation[4] 2, 1,3} ] = {{{0,0,0) . (0, 

0.0}},{{0,0.2),{0,0,0}}.{{0.2,0},{2,2,0}}.{{2,0,0},{0,2,0}},{{0.2,2}, 

{2 f 2,0}} # {{2,2.0}. (2,0 ( 0)},{{2.2,2), {2,0,0}}, {{2,0,2}, {0,2. 0}}, {{O,0, 

1}. {0,0,2}}, {{0,0.3}, {0,0,2}}, {{0,2,1}, {2,2.2}}, {{2,0,1}, {0,2.2}}, {{0, 

2.3}, {2,2,2}}, {{2,2,1}, {2.0,2}>, {{2,2,3}, {2,0,2}} , {{2,0,3}, {0,2,2}}. 

{{0,1, 0}, {1,1,0}}. {{O, 1.2}, {1,1,0} }, {{0,3,0} , {3,3,0}}, {{2. 1.0}, {1,3, 

0}}, {{0,3,2}. {3,3.0}}, {{2,3,0}, {3,1,0}}, {{2,3.2}, {3,1.0}}, {{2.1.2}. {1. 

3-0}}. {{1,0.0}, {0,3,2}}, {{1,0.2}, {0,3,2}}, { {1, 2 , 0 } , {2 . 1, 2} } , { { 3 , 0 . 0} , 

{Q, 1.2}}, {{1,2.2}, {2,1,2}}. {{3,2,0}. {2,3,2}}, {{3, 2,2}, {2,3.2}}. {{3,0, 

2}, {0.1.2}}, {(2,3.1}, {3,1,2}}, {{2,3,3}. {3,1,2}} , {{2,1.1}, {1,3,2}}, {{0, 

3,1}. {3,3,2}}, {{2,1,3}. {1,3,2}}. {{0,1,1}. {1,1,2}}, {{O, 1.3}, {1,1,2}}, 

{{0.3,3}, {3,3,2}}, {{3,3.2}, {3,0,2}}. {{3,3,0}, {3,0,2}}, {£3,1,2}, £1,2, 

2}}, { {*- 3. 2}, {3,2,2}}, {{3,1,0}, {1,2,2}}, {{1,1,2}, {1,0.2}}, {{1,1,0}. {1, 

0,2}}, {{1,3,0}, (3,2,2}}. {{1,3,3}. {3.2,0}}. {{1,3,1}, {3,2,0}}, {{1,1, 

3}, {1,0.0}}, {{3,3.3}, {3,0,0}}, {{1,1,1}, {1,0,0}}, {{3.1.3}, {1, 

2.0}}, {{3,1.1}, {1,2,0}}. {{3.3,1}, {3,0.0}}, {{1,2,1}. {2,1,0}}, 

{{1,2.3}. {2,1,0}}, {{1,0,1}. {0.3,0}}, {{3,3,1}. {2,3,0}}. {{1,0, 

3}. {0,3.0}}, {{3, 0,1}. {0,1,0} }, {{3,0,3}, {0,1,0}}, {{3,2,3}, {2,3,0}} }« 

TupleRaprBsentatian: :uaage - "TuploRepr6flontation[ModuloRing] [Prijnifcive Polynomial] 
designed for fields, e - »- TupleRepirGSBiitatioix[2 ] [ {1, 0 . 0 , 1, 0, 1) ] = ({0,0,0, 
0,1} . {0, 0.0.1,0}, {0,0,1.0,0}, {0,1,0,0,0} . {1,0, O,0. 0}, {0,0.1,0.1}. 
£0. 1,0, 1, 0} . {1,0.1, 0. 0}. {0,1, 1,0. 1} , {1,1, 0,1, 0}, (1,0, 0,0,1} . {0,0, 
1,1.1}, {0,1,1,1,0}. {1,1,1.0,0}, {1,1.1, 0,1}, {1,1.1,1.1}, {1, 
1,0,1,1}, {1,0,0,1, 1}, {0,0, O. 1,1}, {0,0.1.1.0}, {0. 1,1,0, 0}, 
{1,1,0,0.0}, {1,0,1,0,1}, {0,1,1,1,1}, {1,1,1,1,0}, (1,1,0,0, 
1}* (1* 0,1, 1,1}, {0,1, 0,1,1}, {1,0,1.1, 0}, {0,1,0,0,1}. (1,0, 0,1,0} }- 

Minimum Poly: : usage - ^MiiujnumPoly[Moaulolling] [Priml t ivePo lyaaomi al ] [rootPowar] [ 
var] designed fox £ields.e.g, Minimunfoly[2] [ {1, o, 0. 1, 0, 1) ] [3] [at] " 

arae££QMathod: susaga = "Qraef feMathod [GP2minpoly] generates & primitive 

polynoitdal in nmj z4 e.g GrA*f f ertothod £ (1, 0 , 0, 1, 0, 1} ] = {1,0,0,3,2,3} " 

flBquexiceGen^rator; : usage = 
N SeC£ue&ceGexiex^tor[HoauloRing] [Primitive Poly] [mitialCond] [Senl-eASffch] 
e . jj . sequenceGenarator [4 ] ({1. 0, 0, 3, 2. 3}] [{0, 0, 0, 0, 2}] [31] = 
{{2, O, 0. 0, 0, 2, 0, O, 2. 0. 2, 2, 0. 

0, 2, 2, 2, 2, 2, 0, 0, 0, 2. 2, 0, 2, 2, 2, 0, 2, 0}, {0, 0, 0, 0, 2}} 

u 

XnitlalCandifcions: iusage = p initialCondi tiona [PriaitivoPoly] e-o Initialconditiona \ 
(1. 0.0, 3,2.3}] = £{0,0, 0,0,1}, {0,0, 0,0, 3}, {0,0,0.2,1}, {0,0, 0.2, 
3}, (0,0, 2. 0,1}, {0,0, 2, 0,3}, {0,0.2, 2,1}, {0,0, 2, 2, 3}, {0,2,0, 
0,1}, {0,2, 0,0. 3}, £0,2, 0,2,1}, {0,2, 0.2, 3}, {0,2, 2. 0,1}, {0.2, 
2, 0,3}, {0,2,2.2,1}. {0,2,2.2.3}, {2,0,0. 0.1}. {2,0.0.0,3}, (2. 
0,0, 2.1}, {2,0, 0,2. 3}, {2.0,2,0.1}, {2, 0,2, 0.3}, {2, 0,2, 2,1}, 
{2,0,2,2,3}, {2,2,0,0.1}, {2,2,0, O. 3}, {2,2,0,2,1}, {2,2. 0,2, 
3}* {2, 2, 2, 0,1}. {2,2,2,0,3}, {2,2,2,2,1}, {2.2,2,2,3}, {0,0,0,0,2}} « 

GoldSeipidnce: tusaga = "QoldCocjiAanCG [PrimitivePoly] [ luitcoadNujnber] a.g. 
<3olcU5e<xuBxica[{l,0,0,3,2,3}] [33] = {1.0,0,0,0,1,0.0,1,0,1, 
1,0,0.1.1,1,1,1,0, 0,0,1,1,0,1,1,1,0,1,0} - 

Aut ocor r e 1 a t i o»S a gu»nc e : :usd.ge = > *Aufcocoxrxrel&t:ioiiSeiziiance[Gol<i9e^u.eiice] e.g. 
Autocorrelationfiiqueoce [ {1, 0,0.0,0,1,0,0,1,0,1,1,0,0,1,1,1,1,1. 0. 

0. 0.1.1.0. 1.1.1.0.1.0}] = {31.-1. -1,-1,-1.-1,-1,-1,-1,-1, 
-1,-1,-1.-1. -1,-1,-1,-1.-1,-1.-1, -1,-1. -1,-1, -1,-1, -1,-1,-1, -1} " 

Specif ledGoldSaquences: lusage = "SpecifiadGoldSoquepcafltFrinitivBPQly^] [ 
Number-Li ot_] e.g. speci £iadGold5equencafl ( {1, 0, 0, 3 , 2, 3} ] [{1, 2, 3, Last - 

1. Lasfc}] s {£0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,1,1,0,1,1.1. 
0,0. 0,1,1,1,0, 0, 1,1}, {1,0,0,0.0,1,0,0,1,1,1,1,0,1,1.0, 
0.1, 0,1, 1,0, 1,1, 1.0. 0,1, 0,0,1}. (0.1, 0,0, 0,0. 1.0, 0,0,0, 
1.1, 1. 0,0,0,1,0,0,1.0, 0,1,0,1,0,1.1,1,0}, {1,1.1. 1,1,0. 
0,0,1,0. 0.1. 1, 0,0,0,1,1,1,1,1,0,1,0.1.0.1,1,1,1,1}, {1, 
0,0,0,0,1,0,0,1,0, 1, 1,0, 0,1.1,1,1,1,0,0,0,1,1,0,1,1,1,0,1,0}) - 
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AllGpiaSequracess ruwgo = "AllaoldS«qu©ncea[PrimltivoPolyJ e.g. AllaoldSdguences r 
{1,0.0,3,2,3)] s HO, O. O r O, O, O, O, D,0,l,O, O r O, 1,0,1,1, 0,1, 1,1,0, 
O, 0,1,1, 1,0. 0,1,1} , {1,0,0,0,0,1.0,0,1. 1,1,1.0.1.1, 0,0,1.0,1.1. 
0,1,1.1,0,0.1, 0,0,1}, {0.1.0,0.0,0,1,0,0,0,0,1,1,1,0.0,0,1,0,0, 
( i*?'2'^'2'i'2' :i ' :L ' :L ' D) * {!»!' 0,0,0, 1, 1,0, 1,0. 1,0. 1.1, 1.1, 1,0,1, 

V 0,1,0,1,0,0,0,1,0,1,0,0}, {0, 0,1.0,0.1. 0.1.1.1,0.1.1,0.1.0,1,0, 

1.0,0,0.1.1,0,1,0, 0,1,1,1), {1,0,1,0,0,0,0,1,0,1,1,0,1,0,0.1,0, 
1,0, 0,0,0,0,0, 0,0, 1,1, 1,0.1}, {0,1, 1,0, 0,1, 1.1, 1,0, 0,0, 0,0, 1,1. 
2'-*?' 1 ' 0 ' 0 * 1 ' 0 ' L I'LL 0,1.0), {1, 1,1,0,0,0, 1,1,0, O r 1,1, 0,0,0, 
0.1, 0.1. 1.0, 0.0,1,1,0, 0,0, 0,0,0), {0,0,0,1,0,0,1,0,1,0,0,0,1,0. 
1,0,0,0,1,1,0,1,0,1, 0,0,1,1,0,0.1), {1,0,0.1.0.1,1.0,0.0.1,1,1, 
?'2' 1 ' 1 ' X ' 0 ' 1 ' 0 ' 1 * 1 ' 0 '°' :1 * 0 ' 0 * 0 '1'1>. {0,1,0,1,0,0,0.0,1,1,0,1, 
0.0, 1,1, 1.1. 0.0. 0,1. 0,0,1,0,0, 0,1, 0,0), {1,1, 0.1, 0,1, 0,0. 0.1.1. 
0,0, 0,0, 0,0, 0,1, 0,0, 1,1, 1,1. 1,1, i,i,i,o } , {0,0,1, 1.0,1.1.1, 0,0, 

S'J'S'i'^'J-'S'?' 1 ' 0 ^-?' 1 ' 0 ' 1 ' 0 ' 0 ' 1 ' 1 ' 0 ' 3 -)' a. o.i.i* a.0,1, 1,1, 

0,1, 0,0,1,1, O.X.I. 0,0.1,1,0,1, 1,1, 1,0, 1,1,1), {0,1,1,1,0,1,0.1. 




°. 0,0 0,0,0,1,1, 1,0,1,0,0,0.0,1,1,0,1. 1,0, 1,0, o,oJ, 

{0,1 1 01 1 1 0,1,1.1 0.0.1. 0,0, 1,0,0, 1,0,1, 0,0, 0,0,0. 

i' 2' J*^^;^ 1 ;^*?- 1 - 1 '!' 0 ' i' 0 * 0 ' 0,0. 0.1, 1.1, 1,1, 1,0, 0.0,0.1.0, 
S'S' 1 * 0 ' 1 ' 1 *' to.o.i.iri.i. 1,0,0,1, 1. 1.0,0. 1,0,1.1, 1,0.1. 

0,0, 0,0,1, 0.1, 1,0, 0.0, 0.0.1,0}, {0,1, 1,1,1. 1,0. 0.0: 0,1; 

2'i # ?'^i' 0 * 0 *°' :L - 1 ' 0 - o ' 1 ' 1 ' :L ' 0 ' 0 *i' 0 'i>* {1.1.1.1-1,0, 

2' 2* i'?* 0 - 1 '!' 0 ' 0 ' 0 '!' LI. LI. 0.1, 0,1,0. 1.1. 1.1.1), {1. 

O, 0,0, 0,1, 0.0,1, 0,,1, 1,0, 0.1, 1.1, 1,1, 0,0, 0,1, 1,0, 1,1, 1,0, 1,0))- 

C ^i^^^H^^ [Ooldise,] e.g. w. wi fl h to optialoa thie 

™f? ftl vrtM *f fclM** ±fl & bit tranaition . Therefore x (-1) In the x4Bitstat e . it 
aelivex-a tlie two transformed code a . one for each pule© ia either 0 or 1 
co^grranaxoxiDCi.] [{1. o. o. o, o, o, o, o, l, 1*1, o, l, l] ,T, o , £. o. o, o, 1. o, o 
l. i, o. i, o, o, o, o, o, i, i, o, i, o, i, i, i, o, o, o, o, l. o, i i i i o n 

S' 1* ?' S' l* S' S' ?- ?• t ?- i- 1> ?- °- }• *• °- °- i- 2; i; 5; i; S; J; £; 

O, 1. 1, 0, 1, 0, 0, 1, 1, 1. i, o, 1, 1. 0. 1. 1. 1. O. 1. o, o, o, 1, 1, 0. 1 1 o o 

o' a* i' ?- 2* 2- i' ?• °- °- °- x - °- x > x - °' i- 1 i* i i o. i; °: °; 

2' 2' i 0 3" ?' ?* °* °* x ' x * °' °* i- o. i. o. i. o, o, o, i, o, 1. o* i i I i i i 

S' ?' i* i* i- 2' 2' 2' °- °' 1 * °- 1 - °' °' °' °' i- °' S o S; i; i; i; i' J* J- 

0, i, i, o, o, o, i, o, o. o, i, o. o, i, a, i, o. o, i, 1, o. o' ll 1 o 1 1 1 i 2' 

i oioi, o i x o, o o i i, i: i: x - x : x ; 0 ; ?; S; S ^V;^, 1 ' ^ °- 

1, 1, -1,1, I. -1. I, 1, I, 1, 1, 1, -1, —I, 1. I, -1. Z, -X, X, 1, -I. -1, -I 
1, I, 1. I 1. -I 1, -I -X, _x, 1, -I, -1, I, -1, —I, -1, -T, -1, X, -1, I. -1 

I: i i"i "V "t V V I x x ' T - I %- :L v 1 ' a - x - -i^ -i. i, 

I. 1, -I, -1, T, 1, -X, 1, I, -I, —I, 1. ^X, 1, -I, 1, I, -1, -x, -1, -l i i _i 

i. l, ' -i. i« -i. -x. x. -x. -x. x. x. -x. -x. -i, i, —i -d) 
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CodingTr^iBfontlNflW! :\iQaga a "Codingrrrans£oxnNew[X<] [GoldSog] e.g. Wft wlflh to 

optimide this signal wh«a there ifl a sit transition. Tharfiforn x(-l) In the 
xABitState. Xt ddlivers 4 tranaformed codes, one for each pulse is either 0 
or 1 and Chen for the second pulfie CodingTranfl£orm[l<] [{1,0,0,0,0,0,0,0,1. 
1.1,0,1,1,0,0,0,0,0,0,1,0,0.1,1.0,1,0,0,0.0. 0,1,1,0,1. 0,1,1, 1, 
0,0,0,0,1.0.1*1,1,1,0,0.1.0, 0,0, X. 1,1,0. 0,0, 1, 0,1,1,0 ,0.1,0 ,o, 
1,0, 0,1.1, 1,0,1,0,1.1,0.1,1,0,1,0,0,1,1,1,1,0,1, 1,0,1,1,1,0.1. 
0,0,0,1,1,0,1,1,0,0.1,1,1,0,0,1,0. 1,1,0,1,0,1,0,0,1,0, 1.1,1,0, 
1.1,1,1,1,0,1.1,1,0,0,1, 1,0.0. 0.0,1,1,0,0,1,0,1,0, 1.0, 0,0,1,0, 

1.0, 1,1,1,1,1,1,0,0,1,1.1.1.1,0,0,0,0,0,1,0,1,0, 0, 0,0. 1,0,0.0, 
0,1.1.1,1, 0,0,0,1,1,0,0,0,1,0,0,0.1,0,0,1,0,1,0, 0,1,1,0,0,1,1. 

0. 1, 1,1,1,0,1,0,1,0,1,0,1,1,0,0,0. 1,1,1, 1,1,1,1,1,0,1,0.0)] = 
{{X, -1,-1, 1,1,-1. -X, 1, -I. -1,X, -1,1, 1,1, -1,-1, 1,1,-1. 1,-1, -1,-1, 

1, -1,1, -1,-1,1, I, -1,1,1,1,1, I, 1,-1,-1, -I, 1. 1.-1,1, -1,1.1. -I. -1, 

-X, 1, -I, 1, I, -1, X, 1,-1,1, I, -1,X, -1,1, 1,1, -1,1, -1,-1.-1,-*, 1,-1, 
-1,I,-1,I,-1,I.1.I,1,-X,1,-I,1.I.1,-I,-1,I,-1.X. 1,X,1,-I,^1.-X, 
-1,-1,1,1,1, "1,1,-1,-1,-1,1,-1,-1,1,-1,-1, -1,-X,-1,I. -1,1,-1,1, 
-1. -I. -1, -X,-l,X.l. I, 1,-X. -1,1, 1,1, 1,-1,-1, -I, 1,-1,-1, -I, 1,1,-1, 

1. 1. X, -1,1,-1. I. -1.1, -1,-1,1, -I. 1,-1, 1,-1,-1. I, l.-X, -1,-X, 1,-X. 
-1,I,1,-I,1,I,-1,-I,1.-X.1,-I,1,X,-1, -I. -1, -1,1, I, -l.X, 1,-1,-1, 
-I, 1, 1, 1, -I, 1, I. -1, I, -1, -X, 1,-1,1, I, 1, X, 1,1,-1, X. 1,1,-1.1,1,1, 
l,-X,-l,I,-l.I.-l,X,-l,X,-l,I,l,l,-l,-I,-l,l.i. -I, -1,X. 1,-1,1, 
-1,1,1). {-I, -1,X, 1,-1,-1, I. 1.1,-1. -I, -1,-X.l, -X, -1,1,1.-1.-1. -I, 
-1,X, -1,-1.-1, -1,-1,1,1,-1,-1.-1,1. -X.l.-X. 1,1, -1,1, 1, -I, -1, -I, 

-i,-i,i.i,-i,i.i,x,i,-i,-i.-i,i,x,i,-x,-i,-x. -1,-X, 1,-X, -l.-X. 
-l,x,-l,X,l,x,-i.-x,-i, -x,-i,-i,i, -i, i.x, l,i.l,-i,l.x. -1,-1,-1, 

-1,1, -I, 1,1, -1.1, -1,1, 1,-1,1, I. l.X, -1.1,1,1, -1,-1, -1,1,-1,1,-1, 
-I,-1,-I,-1,-X,-1,I,-1,X.-1,-I,1,-X.1,I,-1.-X,1,-X,1,I,-1.I.1,I. 
-1, I, 1, -X. -1,-X, 1,-X. -1,-X, -1,-X, -1,-1, -l.X, l.X, 1,1,1. I, -1, -I. 1, 

x, -i, x, i, i, -l, -x, i, i, i,-x, -1,1, l.x, i.i,i,-x.-i, i, -l.i.i, -x,-i, 

-X. 1,1, -1. X. 1, -X, 1, 1,1, -X, -1, —I,— 1,1,1,1,1, -I, 1, -X, 1, -I, -X, -I, 
1, -X, -1,-1,1, -I, 1,1,-1, -X,-l, -I,-l,-X. -1, -I, -1, -1,1, -I, -1,1, -1, 
-I, 1,1,-1. -1,1,1,1, I, 1,-1}. {1.1.1,1, -1,-1, 1,1, l.X, l.X, -1,-1,1, 
-I, 1,1. -1.-1,-1,1, 1,-X, 1,1. -l.X. 1.1, -1,-X, -1,-X, 1,1, l.X, -l.X, 
-1,1. -1,-1. -l.X, -l.-I. -1, I, -1,1,1, -1,-1, -X, -1, —I, -1,-1,-1, -X. -1, 
I, -1.-1.1,-1,-1, I, 1,-X, -1,1, 1,1, 1.1. -1.1. -1,-X, 1.1, -1,-X, 1,-X, 
-l,I,-l,X,l,X,-l,-I,l.I,-l,I,-l,-I,-l,X,-l,X.-l.-X,l,X,-l,I.l. 
I.1,X,1,-I,-1.-X,1,I,-1,X,-1,I,1,-I,-1,-I,1,-I,1,I,-1,I.1,-X,1, 
X.-1,I,-1,I,1,I,-1,I,-1,-X.-1,-I,1,-I,-1,I,-1,X.-1,I,1,X,1, -I, 
1,-1,1,1,1, -X. -1.1, -1,1, 1, X, 1,-1, -1,-1.-1, -I, 1,X,1,-X.1,-X,-1, 
-I, 1,-1, -l,X.-l,~l,-l, -I,-1,I,-1,X,-1,I,-1,-X,-1,-X,-1,I,1,I,1, 
-I, -1,1, 1,1,1,-1,-1, -I. 1. -X, -1,-1, 1.1, -l.X. 1.1, -1,X, -l.X, -1,1, 
-1,-X, 1,-X. 1,-X. 1,-X, -1,1, 1,-X, -1,-X, 1.-1,-1), {1,-X,1, -X.-l.X.l. 
-X, 1, -X, 1, -X, -1,1. 1,1. 1. -I. —l.X.— 1, -1,1,1,1, —I, -1, —X ,1, —I, —l.X, 
-1,X.1.-X. 1,-1.-1.— X, —1,-1,-1, 1,-1. —I. -l.I,-l,-X,-l,-X,l.I, -1,1, 
-1,1. -1.1, -1.1.-1,-1,-1,1. 1,1, -l.-I.l,I,-l,-X,l,-X,l, -1,-1, -I, 
-l.X, 1, -I, -1,1, l.X, -1, -X, -1, -I, 1,-X, -1.1, 1.-1.-1,-1. -1,1,-1, -X, 
-1, -I, -1,1.1. -I, —1, —X, 1, —X, 1, —I, 1, I. -1, X, 1, —X, —1, —I, —1, -I. 
1.1, -1.1, 1.1, 1,-X, -1,-X. 1,1,1. -1,-1. -I, -1,-X, 1,-1,-1, 
-X, -1.1,-1, I, 1,X, -1,-1,-1, -I. -1,-1,1, -I, l.X, 1. I, 1, -X, 
l.X,-l.rI,-l,-X,l,-X, 1,1,-1.1.-1.1.1, -I, l.X, 1,1, -l.X, 
l.X, -1, -I, -1,1, -1. X, -1, -X, -1,-1,-1, -X, -1, I, -1,1,-1, -X, 
1, -X#l,I#-l # -I # l,-I,l,I,-l.I,l,X.-l.I,l,-i,-l,-X r l.-x. 
1# -X, -X, -1,-1, -I. -1,1.1.1,1.x. 1,1, -1,-X, 1,1, -1,1, 1,1, -1)} - 
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CroaacorrBlationSaquenqa: :usage = " cros acorrela t ionsequence [ {QoldSoql . Qoldseq2 ) ] 
e.g- oeqlist =Bp©cif iedGoldaoquencea [ {1 , 3 , 2 , Q # 1, 1, 1, 1, l f 1* 1} ] [{1,2,3, 
4,5, 6 # 7 # 8 # Last -l.Laot) ] ; Outer [ (Crofl3CQrrBlatiou3equenC«[ {#1, #2 ) ] / /union )&, 
. d^liet^AqliQt^l] = {{{-33,-1,31. 1023), (-65,-33,-1,31,63), {-65,-33,-1. 

f 31,63), {^65, -33,-1,31,63), {-65,-33,-1.31,63), {-65,-33,-1.31, 63) , {-6S.-33, 

^ -1,31,63). (-65,-33,-1,31,63), {-65,-33,-1,31,63), (-33,31)), {{-65,-33,-1, 

31,63)* {-33, -1,31,1023), {-63,-33,-1,31,63), {-65,-33,-1,31, 63 ) , {-65, -33 , 
-1,31, 63), (-65,-33,-1,31. 63), {-65,-33,-1,31, 63) , {-65.-33,-1.31, 63), {-65, 
-33.-1,31,63), (-33,31)), {{-65,-33,-1,31,63), (-65,-33,-1.31,63), {-6S,-33, 
-1,31, 63,1023), {-33,-1,31), (-65,-33,-1*31,63), (-65,-33,-1,31,63), {-65, 
-33,-1, 31, 63), {-65. -33,-1,31, 63), {-65, -33.-1, 31, 63), {-33, 31)), { {-65, -33, 
-1,31, 63), {-65,-33,-1,31, 63) , {-33 , -1, 31) , (-65, -33, -1, 31, 63 , 1023) , {-65, 
-33,-1,31, 63), {-65,-33,-1,31,63), {-65, -33,-1,31.63), (-65,-33,-1.31,63), 
{-65,-33,-1*31,63), (-33,31)), {(-65,-33,-1,31, 63), (-65,-33,-1,31,63), {-65, 
-33,-1,31, 63). {-65, -33 ,-1,31, 63), {-33 , -1, 31, 1023 ) , {-65, -33 , -1, 31, 63 ) , 
{-65. -33.-1,31,63), {-65, -33, -1, 31. 63) , {-65, -33, -1. 31, 63) , {-33,31)), {{-65, 
-33,-1,31,63), {-65,-33,-1,31,63). (-63,-33,-1,31, 63) . { -65, -33 , -1, 31, 63 ) , 
{-65, -33,-1,31,63), (-33, -1,31, 1023), (-65,-33,-1,31,63), (-65,-33,-1,31, 
63) , {-65, -33,-1,31, 63). {-33, 31) ) . {{-65,-33,-1,31, 63) , {-65, -33,-1,31, 63) . 
(-65,-33,-1,31,63) , {-65,-33. -1, 31,63), (-65,-33,-1,31,63), {-65,-3 3,-1,31, 
63), {-65, -33.-1,31, 63. 1023), {-33,-1,31) . (-65,-33,-1, 31. 63), (-33,31)) , 
{{-65,-33,-1,31,63), {-65,-33,-1,31,63), (-65.-33,-1,31,63), (-65, -33,-1, 
31,63), {-65. -33, -1.31, 63), (-65,-33,-1,31,63), (-33,-1,31) , {-65,-33,-1,31, 
63.1023), (-65,-33.-1,31,63) , (-33,31)), ({-65, -33 , -1, 31, 63 ) , {-65 , -33, -1,31. 
€3), {-65,-33,-1,31,63), {-65,-33,-1,31,63), (-65,-33,-1,31,63), 
(-65,-33,-1, 31, 63) , (-65, -33,-1,31. 63), {-65,-33, -1,31, 63), (-33, 
-1,31.1023). {-33,31)), {(-33,31), (-33,31), {-33, 31) , {-33, 31) , {-33, 
31). (-33,31), (-33,31), (-33,31), (-33,31), (-1,1023))) - 



Begin [ - • trivato * ■ ) 



RingDiviaion[ModuloRlng_] [{Denoialnatar., Numerator.) ] z = 
Module [{xl, x2, x3, ae4, x5, x6, Qx, MX), 
Clear [xl] J 

xl [±_, j_] 2 = Hod [i j , ModuloRing] 3 

x2 - Table [xl [Denominator I [l] ] , J] , {j. o, ModuloRlxig: - i)] / 

x3 ^ Complement [ Table [i, {i, 0, ModuloRlng - 1} ] , x2 ] ; 

Qx = {); 

Six c Numerator; 

NoelUpd&te t - 

Module £{uxi, ux2 , ux3, ux4). 

If [MemberQ [x3 , Nz [(!]]], Tnrow[{Qx, Nk. Denominator, False)]]; 
uxl - Take [N*e, Langth[DBnomlnator] } ; 

ixx2 = Position [x2 , uxl[[l]l] //Flatten// <#[[!]] -1) &s 
11x3 = Kod [uxl - ux2 Denominator, HoduloKing] ; 

Qx = tfoin[Qx. {ux2)J; 
ux4 = Drop [Nx, Length [Denominator ] ) ; 
Nx = Join[ne3t[ux3j. ux4 ] ] ; 
Catcb[Tabl«[ Noel Update, (i, 1. Length [Numerator] - Length [Denominator] + l)]j 
{Qx, Nx, Denominator, True) ] 

] 



PropLeadingZerad [Denominator&eq_] ; = 
Module [ (xl) , 
*![<)] {)? 

xl[D3_] /;DS[[1]] = = =0 :~ xl [ Rest [DS] ] ; 
xl [BS_] : = DSs xl [ Denominatorseq] 

] 



JfoBaiblttDivisgra [HoduloRing^] [ordor_] : = 
Module [(xl. x2, x3, x4), 
xl = Table [ (x2 [ij . 0, Modulo Ring -1), {1, 0. order}]; 
x3 = Table[x2[i] f {i, 0, Order)]; 
x4 = Apply [Table, Sequence [ Join [ (x3) , xl] ] ] // Flatten[#. Order] fe // Rest ; 
Map [DrOpLoadlng^Bros, x4] ) 
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PolynomialMultiplication[ModulaRing_] [polyl_# Poly2_] = = 
Module [ {xl. x2 , x3 . x4, x5, x6, xB f x9), 

xl = I,©ngtn [polyl] * iiengt.li [poly 2 ] - X; 
x2 = Tablo[x3*i r {±, °* *1 - 1)3 i 
f x4 = {Tafce[x2, Length [polyl] ] //Rftveree) - polyl j 

xS = (Taxe[x2, length [poly 2 ] ] // Revarne) • poly2; 

x7 / Tablo^oef£±c;±ont[xfi. ac3 * (xl - i - 1) ] /V Mod[#. ModuloRing] (i, 0. xl - 2) ] r 
X B = tfoin[x7, (»od[x6 /. x3 -> O. ModuJLoRinH] } } // DropLBadingZerod] 

Wc define polynomial muhplicaiion, modulo a primitive polynomial 

ModuloMultplication[Haduloaiiig_] [primitivePolynomialJ [polyl_, poly2_] : = 
Module [ {xl, x2 f x3, 3C4}« 

xl = polynomialMultiplication[HoduloRlng] [polyl, poly2 ] ; 
ac2 = RinflDiviaion[Moeiul6Rlng] [ (PrimitivePoJ-ynomial, xl)]; 
x2[[2J] // DropLaadingZeroa] 

RlngPowear [RiugModulo_] [PriinitiveElwiantJ [Element_] s = 
Module [ {jcI, et, NoelUpdata) # 
xl ~ Length [PrinitiveBlttment ] -1; 

fit = {i>; 
Noel Update = = 

st = MaduloKUitplication[RlxigKodulo] [pT*lmitivaElementl [ at, Element]; 
Tahle [Vodlupdata, {i, 1, 2* 1 -!)]] 

CycllcMul t iplat lveGroup [ModuloRlng_] [P:rlinitivaPolyrvoin±al_] = = 
»odule[(xl, x2, x3. x4, x5, x6 # x7 # x8, x9), 
xl =^ (t,engtn[PriJnitivePolynomial] +1) /2 //Floor; 

X2 = FOflsiblaDiviSQrfl[HoduloRing] [xl] ; 
x3 = MaptRingDivlsion[ModuloRincr] [{#, Primi tivttPolyuomialJ ] 4, x2] ; 
= select [x3 , # [ [4] ] == Tma &] ; 
xB = $Qlect[x4, #[[2]] == (1) t] //TraneposB// Join[#[(l3], #[[3)]]tj 
x6 = Selact[x4, #[[2]] == {ModuloRing - 1} fc] // Transpose // Join[#[ [1] ] # #[[3]]]&i 
x7 = Join[{{!}, (ModuloRing -1))* *5, x6] // Union; 

xQ = Map [ModuloMUltplicationfModuloRing] [P?:imi t ivePolynomial] [# # *]*, x7] // Union; 
x9 a RingPowor [ModuloRing] [ PrlmiCivoPolynomial] [x6 t [2] ] ] // RotataRigfct ] 



ZaroSequances [ModuloRingJ [order__] - = 
Module [ {xl , X2, x3 r x4}« 
Xl ^ Tabl«[{x2[i] f 0, 1>, {i, 1, Order}]; 
x3 = Table [ ModuloRing / 2 x2 [i] , {i. 1, Ordar }]; 
x4 = Apply [Table, seguenco [Join [ {x3 ) , xl] ] ] // Flatten [tt, order - l] & ] 

ZerQ?ad[order_] [Element^] s = 
Module [ {xl, x2, x3}, 
xl = Table [0, {i, 1. Order - Length [Element) ) ] ; 
x2 ^ Join[xl f Element]] 

TwoAdicExpanaion [HoduloRing_] [PrlnitivePolynomialJ :- 
Module [{xl, x2, x3, x4, x5, x6, x7, xfl, x9) , 

xl = LsngCh[PrimicivePc»lynoaial] - 1/ 
x2 = cyclicMultiplativeQroup[ModuloRi»g] [Primitive Polynomial] ; 
x3 = Hap[ZeroPad[xl], x2] ; 
x4 - Join [ {Table [0, (i„ 1, xl}]}„ x3] ; 

x5 = Function [x, Hod[x[ [1] ] + 2x[[2]] r ModuloRing] ] ; 
x6 = Table[{x4[[i]] r x4 [[j] ] , x5[{x« ( [1] ] , x« [ fj] ] }] // DropLeadingZeroe], 
{i, l r Langth^]}, {j, 1. Length [x4] } ] // 
Flafcten[#, 1] &] 

a [HoduloRing_] [Primitiv^Polynojnial^] [TwadicElemantJ ; = 
Module [{xl, x2 # x3, x4, x5, x6, x7>, 
{xl , x2. x3) ~ Two ad ic Element j 
x4 = I*ensrtn[ PrimitivePolynomial] - 1; 

x5 ~ ModuloMultplication[UoduloRing] [PrlmltivaPolynomial] [xl, xl] // ZaroPad[x4]; 
x6 = ModuloMultplic at ion [ModuloRing] [PrlmitivePolynomlal] [x2 , x2) // Zeropad[yA]; 
Mod [2 xfi -r x5, ModuloRing] ) 
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AutomorptiismSigma [Modul©Rixier_J [PrlnitivePolyxumiialJ r = 
Mpdulfi [ {xl, x2, ac3, x4, x5. x6 f x7 # x8, x9), 

xl = TwoAdicBxpaiieicm[ModuloFiing] [Primitive Polynomial] ; 
x2 = Map [cf [ModuloRing] [Primitive Polynomial] , xl] ; 

x3 = Map[ZaroPad[&ensrth[Pr±mitiveFolynoxiiial] - 1] , Tranopoflft[xl] [ [3] ] ] ; 
{x3 , x2) // Txandpoflo] 

UnitflRing[ModuloRing.] [PrimitivoPolynomialJ = = 
Modula[{xl f x2, x3. x4} # 

xl = A.utoxnorpbJ.cmJSigma. [Modulo^ing] [Primit ivePolynaml&l] i 
x2 - xl // Transpose // #[[l]]St; 
x3 = X2 // Mod [2 #, 4]& // Onion; 
Ccmplemant [x2 , *3 ] ] 

UniteR±ng[41 [{1, 0. 0, 3, 2, 3}]; 

r[Ordor_] [ <r ] [x_] t = (KfBatLiflt[a, x, Oxrder - 2 ] // Apply (Flue, #]&) 

iTrac eRepresentat Ion [ModuloRing_] [PjTimitivePolynomial__] i = 
Module [ {xl t x2, x3, x4. xSj x6, x7 # x8, xS> } , 
xl = Longtb [PrimitivePolynonLlal] - 1; 

x2 =i AAitomorpbigmBigma [ModuloHlng] [PiriiiiitlvePolynomial] ; 
Map [ (x3 [Bvaluato [tf [ [1] ] 1 ] : = Evaluate [# [ [2 ]]]) x2]; 
x4 = Map[{#, T[xl] [x3] [#]}&, Transpose [x2] [[!]]] //Modfft, woduloRlnglt 



TupldRepresentation [KoduloRiAg n ] [Frimit ivoFolynomi&l_] ;= Module [ (xl, sc2, x3}, 
xl = Length[ Primit ivoPolynomial] ; 

x2 - NedtLlst [ModuloMultplication^ModuloIling ] [PriniitivePolyxuimiall [ (1, 0), 

{1>. 2* (xl - 1) -2]; 
X3 = Map[Z«roPad[xl -1]. x2]J 

This algorithm docs noi always yield the correct value for the Minimum polynomial. 

MlntmuaPoiy [ModuloRing_] [ Primit ixr€»Pcilynom4.al_] [rootPover_] [var_] s = 

Modul8[{xl, x2, x3 # x4 f x5, x6 , x7 # xB, x9, xlO, xll, x!2. x!3, xl4 # xl5, torn} # 
xl ~ Length [ Pr i mi tiirePolynoiiiial) -1; 
x2 » 2* 1 - 1; 
x3 = Moetliist [Mod[# + #, x2]& # root Power, x2] ; 
x4 = Position [x3, rootPowor] // (# [ [2, i] ] - 1) t; 
xti s Take [x3 , xd ] s 
x7 = T\iplBRopresenta.tion[ModuloRirig] [PrimitivePolynomial] ; 
X0 = Tablafi, {i, 0, x2 - 1} ] ; 
x9 = {xQ, x7} // Transpose; 
xlO = Map[(var + tom[#3)&, x6] // Apply [Times , #]&//£xpand; 
rule = tom[n_] tom[mJ -> torn [n + jn] ; 
xll = xlO //„ rule? 
xl2 = Table [ Coefficient [ xll, var, xl - i], {i, 0. xl}]; 
Map [ (torn [Evaluat6[#[[l] ] ] ] = = Evaluate [# [ [2 ] ] ] ) & , x9] ? 
com [n_] i = evaluate [ torn [Hod [n, x2 ] ] ] ; Mod [ 
Xl2, HaduloRlng] ] 

Gras f £ dHathod [GF2minpoly_] t = 
Module [ 

{xl, x2odd, X2avan # x3odd # x3aven, x4aven, x^odd. xS # xGeven, x6edd, x7 , xB , x9 > , 
xl - Length [QF2min©oly] ; 

x2odd = Table [{1. 0}, {i. 1, (xl +1) /2 //Floor}] // Flatten // TaJca [#„ -xl]i; 

x2even = Table [ {0, 1), (± # i, (xl + l) /2 //Floor}] // Flatten // Take[# ( -xl]fc; 
x3evan = GF2ninpoly x2even // DropLd&dingzeros; 

x3odd = GF2raixipoly x2odd // DropLeadingZerosi 
x4even = PolynomlalMultiplicatloii[4] [x3even # x3even] ; 
xAodd = Polynomial ami tipl lea Llon[ 4 ] [x3odd, x3odd] ; 
X5 = Max [ (Langtb [x4even] , Langth [x4odd] } ] ; 
xfieven = ZeroPad[x5] [x4even] ; 

x6odd = ZeroFadfxS] [x4odd] j 
x7 = xfievan - x6odd; 
xQ = (x7 Sign [x7 [ [1] ] J ) //ModC*, 4]Bk; 

Join[x8 4 {0}] // Partition 2] // Vranepaee [#] [ [1] 3 & ] 
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S^enc^neratortKoduloRingJ [PrimltivPolyJ [IxxitialCond^] [Se«L«mSFfc*0 = - 
MoflulB[(xl > x2Stato, x3Update, x4>, 

xl ~ Mod [-Kept [PrimitivePoly] , ModuloRing] ; 
x2State =* ixw-tlaXCoxxd. ; 

x3 Update == Hodul«[{), x4 = Laat [3f23tate] ; in . vA1 . 

x2£tafce - Join [{£*<* [estate . xl. ModuloRins] > , Drop[x2State, -1]]; x4],- 
{T&ble[x3Update # {i, 1, SaqLoiisrtfa) ] , x2st*fce) ] 

InittalCondlt iona [PrlnitivePoly_] = = Module [<*1, x2, x3, x4, xSJ, x6, x7), 
xl = Length [PrlmitivePoly] -1; • . Min . 

x2 = «abl»[<x3[il. D • 1} # {±. 1* »U 3 * = * ablo r [ J f 3 ^ 1# }1 ' 

x5 J = (Table @0 Sequence [ Join l{x4>. x2]]) mtt©n[tt, xl - X] &; 
x6 = Mod[x5J[ [1] ] + 2xSJ[[2]]. 4]; 
x7 = Map[ Mod [xStft [2]] + 2 #, 4]** x5 JJ ; 
Join [x7 , {x€}]3 

GoXAfleciueiica [PrimitivePoly^.] [IiiifcCoiidJaunibex-^j 7= HoduleK^cl, «2, x3, x4}, 
xl = Length [pritnitivePoly] -1; 

^2 ^ initialConditioMlPrimitlvePoly] [ [XnitCortdWuinber] ] J 
x3 - saquBnceGenerator[A] [PrimitivoPoly] [x2] [2*xX - X]; 
x4 = x3 [ [X] ] /. {0 0, X 0, 2 X, 3 1) ] 

AutocorrolBtio»36quertce[GoldSe(iiiexice_] = = 
Vddule[{xl, x2 State r x3US»dato, x4}, 
xl = GoldS sequence / . 0 -» - X j 

x3UPdate"=* toduXeCO, »4 = xl . x2stat«; Estate - RotateKight l x2Stat B ] ; xA]; 
CPalaXe[x3 Update, {1, X. Lomotil [ Go XdSenuexice] >] 1 

Speci£±edQoldSa*iueiiGoo[PrlTOit±ve£'oly_J [NunbartistJ = = 

Module [{xl, x2, ac3, x4), 
xX = tengtliCPriwiitlvePoXy] - X* 
ac2 - BJumbertiofc / - Last -* 2 A xX ♦ Xj 

Hap [GoldEequancQ [FximitivePoly] *2] 3 

jaiOoldSoqueAcea [PrimiCivePolyJ = = 

ModuXet {xl f x2, x3, x4), 
xl = LexigtH [primitive Poly] - X; 
X 2 » TaM.e[±. £i. 1, 2*xl + 1)]? 

Map [GoXaSa<*uenc;a [primitivdPoXy] [*]&• x2] ] 

CroascorralationSeiiaence [ {Goldfieql^y GoldSen2_} ] : = 
Module [{xX, x2stata, x3npdat B , x4), 
xX - GoldSeql / . 0 -X; 
x2State = ( GoldSeaZ / - 0 -» -1)? ^ r _ _ ^ _ . 

x3 update ModuletO, x4 = xl - estate; x23tata = RotataRlgbt [x2Stat6] ; x4] , 
Table [x3 Update, (1, X, Lengtll [Cold sea; 1] >] ] 

CoOingTranafona[L_] [GoldS*q_] s= . 

Modulo [ {xl, x2. x36tat a< x4Bitstate, xSAccfltate, x6Update, cl, c2, xvj, 

x2 = Goldseq / . 0 -> -X; 
x36tate = -X; 

x4BitSbate = (-1) * liaefc[x2J; 

x5Accstata =0; „ L 
x6Update [x_] s= Module [{), xSAcCBtato = x ♦ xSAccStata; 

x3Stato = x3State * (-1); cl - x3stat«I A xSAccState; 
c2 = x36tate I A (x5AccState - x4Bitstats) ; 

x4Bit£tato = x> {cl, e2} ] ; 
x7 = Map[x6Update, x2] / J Traadpose ] ; 

CodixxgTr-axisforTiBJew[i*_] [Golds : a 

Module [{xl, x2, xSState, x4BitState, xSAccBtate, x 6 Up date, cl, c2, X7>, 
x3 - OoldSeq /. 0 -> -1/ x3 = CodingTranB [L] [x2]; 
x4 c Coding^rans[L] C-x2] 5 
{X3 [[!]], x4 [ [X] ] , x3 [ [2] ] , x4 [ [2 J ] > 

] 

CodingTrana[L_) [siPolarBit3*g_] = = Module [ {xl, x2, x3state, x4BltStata. 
xSAccState, x6 Update, cl, c2, x7 > . x2 - &iPoXazrBltSe<(; x3 State = -X; 
y4BitState - -t.astC^2J; xSACcState =0/ x6Updata[x_] t= Module [ { > , 

xSAccstata = x ♦ xSAccStatfl; x3State = x3State*-l; cl = x3 State I *x5AccStatej 
c2 = x33tate * I* (xSAccState - x4BitState) ; xABitstato s= x; {cl, e2)]; 
x7 = Tranepnae [xfiUpdato /§y2]] 



P 50 R-780 Job-222 

FAX:00 44 1276 677720 P. 050 

10 



11-03-98 16:41 00 44 1276 677720 

l-Mj^98(WED) 16:53 NMP PATENTS UK 

Rin s FunctionsDiffEncodedParenT2. nb 



( 



End[J 



11-03-98 16:41 00 44 127S 677720 p. 51 R-780 Job-222 

MAR. ' 98 (WED) 16:53 NMP PATENTS UK FAX: 00 44 1276 677720 ] 

LaurcntFunctiansPatenr2.nb ( r 



1 



r 



BeginPackage [ ° Laur datPvinc t iona * w ] 
Ts :usaga = "This is the symbol period- 
BT: : usage = "This ie the usual product" 
h;:usage = "This Is the raw gaussian pulse" 
9::usage = "This denoted modulation index Pl M 

:ua&ge = **iJ/[L, t] is Uurento function" 
hPiltBred: :usag6 = -This Is the filtered gauflaiaA pulse" 

Pha B eAn 9 lB:tu B age=-Thi a function taKes some time t° e^ulat. . *ta* . following 

code will drav a graph, of the £unct3.onPhaaoFuncticMi[t_] -M L <Pi.,c J r 
phasepoints =Table [ {t,PhaBaPanction[t T] } /^N, {t , 0,1- ,1/40}]; 
LifltPXot [pliaflepoiiftCB, Plot Joined -> True]" 

FhaedAngi«Fa S t: :u B age -PhaBaAngleFaettL] [t] speeds up the calculation 
ofThafieAngle ^calculating Phae ©Angle [1.] [tj with X, numeric a~d t 
symbolic. It takes some time to calculate ► ■ 

£::U8age= M S = Sin[*]° 

j*;usage- "J = fa* " " 
C:iuBage= "C s= cofl[»] p 

Mr tupage ~ "M = 2 t - 1 " 

Modulationlndext : usage "Modulationindex = h» 

LaurentBi :uflage = -L B uir«nt3(l.] [n] [t] - Sin[ *[t.,t + n T] ] / S" 
LaurentC = : usage = -LaurentC[L] [K] [t] 13 Laurente C*, t - 
AlphaKl: ;usage = 'AlphaKl [U<] [K, i] is taurente ct*..i " 
LaurentLK::U8aga = H t.aurent:LK[L] [K] tfivee the support o£ 
The stair of Modulator Definitions 

A ^^"tia"t a ^ sets u * 

of parlor states of A*,* that the modulator went thgough to get to the 
conetelation point specified toy AccumulatedPhaae which j-b really A 0>0 - 

AKH: :ueage = -akn[l] [KJ [ {State, AccumulatedPhaae} ] defines in terms of A Q ,„- 

Kodul at ingPul aft: :uaage = , 
-The puase ie assumed to have the following structure Pulse [&] [K] [t] ■ 

NUnberOfcurvBO! :usage = "The number of pulses used by the modulator" 

S ^p?in^uCri;i"2^the interval between ea.pl,, of the output of the modulator- 

InitialState: luaage = -InitialStato is the set of hits that are assumed 
to be present bafora i.e { , a.j , » - - > " 

StartingQuadrant; :uaage = » st ar tingQuadrant = A G . -i ana is a number" 
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Modulator: :Ufl&ge ~ "Modulator [Ii] [BitSaq.Opta] asBumad the following 
default: options StartlngQuadrant -» 0 f InitialState -* Table[l, {x r 
1 # 20} ] , Samplicgliitarval -* T/32,MumbarOfCurve0 -» 4 ,Wodul a tin g^uls e^ 
^ Laurence - The Pulse Is assumed to have the following ptructuro Pulse [!■] [K] [t] ' 

Start of the Receiver Functions 

FiltPulflAMuaage » "Tha default pules and Is called by FiltPulse [I*] [X] [t] " 

SyncS^nple: tusage = "Given th&t the sajupling interval is tben eync 

sample haa tang -16 to 16 and this moves the point ut*ed to demodulate " 

Receiver: tusage = ** Receiver [Ii] [Xnput Sec;, Opts] asaumes the following 
default options { 3tarting0uadxant-»0 ^ Initial States {1, 1, 1, 1, 1, 1, 1, 1, 
1, 1,1,1,1,1, 1,1,1,1, l # l} f Samplin£finterval -> T/ 3 2, Modulating Pulse — 
FiltPuiae, Sync sample -> 0 , NunberofcurveB -> 2 

Receiver Propers z us age = "ReceiverProper [L] [ StartingQuadrant, Inibialstate, 
SamplingIntarval f ModulatingPulBe, SyncSanvle^NuiabQrOfCurveo, InputSaq] 
ie called by Receiver af bar all the options have been resolved" 

Begin l Private ] 
VtO0[2) 



h[t_] == 



2 jtBT 



V 2 tt a T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However, this version of 
Mathemaiica gives an error. 



T 

h< lx [t_] = = Relea»a[Module[{r} f J"^ ^ dlt] ] 



~r h[fc - t:] 

~r 

h7iltorad [Pulfi«Wldth_] [t_J s = Module [ {xli X2, x3} , atl = 

* r^, . r_ . r r- h[tl~r] , , Pulo*Width PulseWidth T % ,, 
(N[#l. 40]&) [Tahla[{tl, J ^ — i— dlr}, {tl, , m — }]]? 

x2 = laterpolation{xl] ; 3c2 [t] ] 

9 : = n [ModulaelonzndexTT] 

c t = cos[tf] ; 
£ : = Sin[ff]; 
J == ( ■ // chop) i 
M ;= a*"*, 

PhaaeAngle [t_] /; t £ 0 == 0 

PhaaeAngle [t-_] [t_] y ; tiLT : = s 

PhaseAngle [£■_] [ t_] : = 

Phaa©Angle[L] [t_] = Module [ {xl, x2, k3, x4, x5, x6} , xl =Wiltarad[3LT] [ tl - L T . ] - 

2 1 

x2 =Table[{t2, * | Evaluate [xl] dtl}, {t2, 0, I* T, — ^— }]; Interpolation [x2] [t]l 

J.LT 100 ^ 

We need to ensure that we always cumulate PhascAngle with t symbolic first lo do the calculation only once 

PhaHeAnfllePaat[LJ [t_] : = Modul©[{xl), xl[tt_] = PhaaeAn$la[L] [tt] ; xlffc] ] ; 
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t — ] /: Q < t < tiT i= FaaseAn^lBFast [I,] [t] 
tML_, t_] /; 2liT > t * LT : = D - PhapaAnglfiFiflt [L] [t ~ LT] 

iieed to put this to avoid the function being extrapolated 

r t_] /; 1 <0 < t < LT) &fi I (jLT>tiLT) :=0 
LaurentS [u_J [t_] : = Sin[ ^ t + n T] ] / S 

AlpnaKX [K_, i_] /; (0 < i < && (o <= k < 2"- x ) - = 

Modulo I {xl, x2, x3, KNUM) , xl :a |x2 = Mod[KNuni, 2] ; K*Tum = — — ; x2j ; 

Kttum = K; x3 = Table [al, (i±, 0, I.Z. - 1) ] ; x3 I±J ] 

LaurenbLK[Ii_) : = 

HodulB[{xl} , xi m Tafc>lo[i. (2 - AlpnaKi[ii] [K, ii] ) - ii, {ii, 1, I, - 1} ] ; Mia [xl] J 

touureutC ft_] /; O < = K < 2** :^ 

Liaure»tS[I«) [O] [tj LaurantS [L] [ii + t JLlpH&KX [Z.] [X, ii] ] [t] 

11=1 

The start of the modulator function 



ANKInitialStateSetTTp [K_] [InitBitSBQ jf AccmaulatedPhaseJ : = 

Module [ (xl, x2, x3 # x4, *5, actiphaaft. initbitse<3) , 
Inifchitootx = XnitBitSeq; 
acupnase = AceuinulatedPbaso; 
T^idateSeq s = 

Modulo [ { ) , xl = acuphaa^ - Sum [inithitaeq [ [i ] ] AlphaKI[L] [X, i] J {i, X, r. -!>]? 
acuphaaa - acmphaaa - Firat [initbitfleq] ; initbltaea = Rest [ initbltseq] : «1] ; 
Table [ Updateseq, {i, l, Laurenti.K [I.] [K] ) ] ] 

AKM[L_] [K_] [{3tato_ f AccumulatedPhaae J ] : = 
AecumulatedPhase - Sum [State [ [i + 1] ] AlphaKI[L] [K, ±] , {i, 1, L -i)J 

OpbiOnd [Modulator] := {StartingQuadrant -» O, initlalState -* Table[l, {i, 1. 20)], 
sainpllngmtorval T / 32, Nuinborof curves 4 , Modulat^gPulge -► Laureate} 

Madulator[L_] (Bitaecr * Opt a ] r = 

Module [ 

{xl # ac2. k3. ac4, acS, xfi, fltate, AcciunalatedFhaae, a*q, AJKHState, Curves, Pulse) r 
xl = sanplinginterval / . {Opts} /. Options [Modulator] j 
state - InitialState /. {Opts} /. Options [Modulator]; 
»3 = StartingOuadrant / . {Opts} /. Options [Modulator] ; 
X4 = Sanplinginterval / . (Opts) / . Optiona [Modulator] ; 
Fulfi* = ModulatingPulaa / . (Opts) / . Options [Modulator] ; 
Curves = (MumbarOf Curvaa / . {Opts} / . Optiona [Hodulatdr] ) - 1; 
AccunmlatedPhase = x3 ; 
seq = SitSea; 
Tabla [ 

AKNStata[K] = ANKInitialStatofletDpfL] [K] [atata. ACCumul&tadPh&aA] , {X, O, Curvas) ) ; 
x5 : = Modulo[{>, state = Join [ (First [aeq] ) , oroptatata^ -1]]* 
AqcumulatedPaase - AccunulatedPhaso ♦ First [soq] ; 
seq = Rest [aaa]; 

Table [ AJOffSt ate [K] = Joinf [AKN[L] [K] [ {state, AccumulatadFhaao) ] > . 

DroplAXMState[K], -!]]# {K, 0, Curves)] j 
x6[r_] = sum [ sum [ (J) »»««•[«] [ [t - in p„ lfiA[t-] [ K ] [t + ±t), 

{i. 0, LaurontLK[L] [K] - 1) ] , (K. O. Curves) ) ; 

Table[x6[t], {r, 0, T - *4)]]; 
Tatole[x5, {Jcfc „ 1, Length [Bit eaq] ) 3 // Klatten] 

Options [Receiver] : = {StarcingQuadrant 0, 

initialstate - {1. 1, l r 1, i # l, i, i, 1. 1, X, l, l. l, l, 1, 1, 1, i), samplinginter 
Modulat ingPul a e -* FiltPulfla, Gy&cSample -> 0, NiunbarOfCurveB -> 2); 
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Receiver [inputSeq.. Optu 3 ; = 

Module [[xl,x2* x3, x4, x5, jc6] f 

xl = StartingQuadraxLt / - {Optfl} /. Options [Receiver] ; 
x2 = initialStata / . (Opts) / . Options [RecaivBr] ? 
x3 = Samplinglxiterval / . {opto} /. Options (Recaivor] ; 
x4 = ModulatingPultie / . {Opts) /. Options [Receiver] ; 
x5 = SyncSanpla / . {Opts) / - Options [Receiver] ; 
x6 = NuraherOfCurvefl / . {Opts} /. Options [Receiver ] 3 
Receiver Proper [I,] [xl, ac2, x3 , x4, x5, x6. Input Beq] ] 

RecaivarProper[L_] [ StartingQuadrant.. in±tial£tate_, Samplinglntflrval., 
Modulatingtulsq^, SyncSampl*.., NumberOf Curve s_. Inputs eg ] 1 = 

Moaule[{xl, agn, RoceivedSeq, ExpectedValue, aecf. Receivettext, 2D, 
xl - Ty 9amplinglttl:erval; 

ReceivedSeq = Partition [Inputs eg> xl] // Tranapofla // # [ [syncsample + 1]]&> 
Expaetadvalue = 

ModulatiagPulse J [O] [ (LaurantLKfL] [O] / 2) T + SyncSampla Sunplinglntarviil] ; 
& - StartingQuadzaat; 
agn = O ; 

ID = {); 

flaq ~ Racoivadsaq; 

ReceivaRaxt : = 
Module [{xl, x2), 

= ( {-l) -ffB J* Firat[ae«]) // Im; 
Xf[Aba[acl - ExpactedValue ] <= aba [xl + ExpectedVAluo ] - 

r> = Join[ 3D, {1}]/ A s JT + 1, 3D = Join[ 3D, {-1>]j J! = J? - 1 ] ; 
sea = Rest [eeg:] ; sgn = Mod[egn 4- 1, 2]]; 
Tanle [ReceiveNext. (i, 1, Length [ ReceivedSeij] }] ; 



ID] 



Ewd[] 



EndPacJcA^a [] 
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Meeds [ * WtiXitiaa * Notation* » ] 
NotebookOp«n [ " LaurontFale t to . nb r ] ; 
Not ebo okopen [ " Laur antPalette2.nt" ] ; 

Symbolize [a] 

NotatiantL^,^ « L&ur©ntLK(L_] ] 

Pa ctern : :0at3ym - First element in pattern Pat tern [ B , _] ia not & symbol. 
Syrot>olia:e [bf n ] 

Notation [h £lx [*uaaeWidtii_] Jayilt«red[I»ulflewldtiO [fc__] ] 

Notation « PliaflaAnoiB [t_j ] 

Notation[o Ll >Kj> ii <=> AlptiaXX [EiXi_] [K_, iij ] 
Needa [ "Calculus * L<ap3.a.c^T:r«Tifi fonn x m ] 
Needfl [ -bosael * - ] 
Needs ["Convdlve * " ] 



11-03-98 16:41 00 44 1276 677720 p. 56 R-780 Job-222 

Mjf*98(WED) 16 = 54 NMP PATENTS UK FAX = 00 44 1276 677720 



besseLm 



C 



Needs I "Calculus % LaplaceTranaform* ■ 3 f 

BosbqI * Private* InvftraetaplacoTranaform^InvBrfloLaplaceTranflformj 
BegixiPackage [ "bessel * " ] 

BasselFilter: :usage ■» 

"BdsselFilter [order] [LowPaaa [Cutoff Frequency (AngularFrequency) 3 ] [FrequencyResponse] [a] 

generatfte a filter witli a. 6dB point specified.. To get a 3dB epecificatian 

use 

Be^flalFiltar [order] [LowPass3dfi t Cutoff Frequency (AngularFrequoney) 3 3 [FraquencyResponse] [ 
s 3 - * • • 

To get firont s plane to the £ domain make the transformation s — > 2 Pi f I " 

GrDupDelay; :uflage = H GroupD©lay[Croq] [e] [Angle] gives the group delay seconds at 

the specified frequency . » . N 
IiowPaad usage « ■ Lowpaad - . " 

LawPaB83dB: :uflage = "To get a 3dB point: far- the sesaal filter " 
FrequencyResponse: :usage M Continuous 11 
inpuleaRaflponee: :usage ^ ■* impulse Reap" 

DigitalFreguencyReaponao : :usage ""Digital Frequency Response " 

FrequencyRoflponeePunction: ;uaafle= p Functional form of Frequency Response function" 
RequeptadFrequancy: :usage b "Used in GfOUpDeXay" 
Begin [ » * private » ■ 3 

d[n_] [fc_] <2 n - k) 1/ {2* (n - k> (k! ) (n -k) l ) 

Bessel [nj [e_] : Sum[d[n] [fcj a*]*, <k,0,zi>3 

Definition of 1-ow Pass Filter ■ *) 

BeoselFiiter [order_3 [i*owPass3dB [cutoff Frequency_] 3 [FraquencyRaflponse] [o_] z* 
Module I ( e # t «npf il tar , xl, x2 , x3 f Omega , f req3 dB , ord) , 
tempf ilter [ordL_] [epj := d[ord] [03 /Bessel [ord] [ss] ; 

amplitude = tempfilter [order] tl Omega 3 tempf ilter [order 3 [-1 omega] //Simplify , 

£req3d» «= Select [Omapa / . N [Solve [amplitude l/2,omaga] ] # Ai>s C#l == ft&] [ [1] 3 ; 
tempfilter[ordar] [s/Cutof fPraquency £req3dB3 1 ; 
( *BeaeelFiltar [ordar M ] [IiOwPasa3dB [Cutof f Frequency.,] 3 [FrequancyRasponee] [a_3 1 = 
Module [{as, tempf iltar,xl, 3t3 , x3 , Omega, £req3dB, ord] - 
temp£iltar[ord_] [ea„] - d[Drd] [0] /Beeaal [ord] [bb) ; 
amplitude - Simplify [ComploxExpand [tempf liter [order 3 [I 
Omega 3 Conjugate [tempf ilter [order] [T Omega] 3 3 3 ; 
x3 a Numerator [amplitude] ; 
xl t= Denominator tamplibuda] ; 

x2 = Simplify [CompleXExpand [xl/ . Power [Alia [x_] , 4] -> Power [Re [x] *2 + 
im[x] A 3,2] 3 1 7 

amplitude - Simplify [x3/x2 3 ; 

freq3dB = select [omega /. N[ Solve [amplitude == 1/2 „ Onega] 3 „ Abo [#3 == #t] [ [13 3 j 

temp filter [order] [fl/cutof fPraquency £req3dB3 3;*) 

BeaaelFilter[order_3 [LowPaflfl3dB [Cutof fFrequency ] ] [impuleoRe spouse] [t 3 :« 
Hodule[{B r tt # kl} ( 

xl = N [BaaaalFilter [order] [LQwPaaa3dB [Cutof fPraquency] 3 CFrequencyResponse] fa] ] s 
InvertieLaplaceTranaf orm [xl f a r 1 1 ] ; 

BasselFiltar [order_] [LowPaaa [Cutoff Frequency.] ] [Frequency Response] [a ] : t= 
Block [ (ss, tempf ilter >„ — 
tempfilter = d[order][0]/ Apply [Times, 

(a - as/ [d [order] [0] * (1/ order) ) ) /. 
N [Solve [Bessel [order] [aei] == 0, sa]]]; 
tempf ilter/. s -> a/CutoffFrequency 
]/; 0<= order < ,15 

< * Definition of BandPasa Filter *) 



BeaselFllter [order_J [BandPaes [Lower_ ( Upper_] ] [FrequencyResponaa] [s_] : = 
Block [ { as, tempf ilter] , 

tempfllter » d T order] [□]/ Apply [Times, 

(s - aa/ (dtorder] [0] * (1/ordar) ) ) /. 

N [Solve [Bessel [order] [ae] == 0, ss] 3 3 ; 

tempf liter / , 0 -> <a*2 + Lower Upper) /(a (Upper - Lower)) 
3/; 0<= order < 15 

BaaaalFilter [order J [LowPae m [Cutof fFrequency.] 3 [ImpulaeResponse] [t_3 :■ 
Block [ (83| tempf ilter* denomf ilter, numfiltar, roots) , 

tempf ilter = (c - as/ (dtorder] [0] « (1/nrdar) ) ) /. 

W [Solve [BeacJel [order] [ea] 0, aa] ] ; 

tampfUter = tempf ilter/ . s ~> a /CutoffFrequaney; 
tempf ilter » Simplify [tempf ilter] ; 
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r Title : *) 
(* :AUtbora: *) 

J r context : Global » *) 

(* : History: Started on IS «mne 91 
Worlied on* > 

(* z Functions = NoeJ-Convolve *) 

(*■ : Sovlztco : *) 

( * :Waraing: *) 

( * ;Hath«ABtica Vftraion: 1.2 * ) 

(* =Ii imitation: * ) 

C 

:DiBCuadio&! 

the program calculates rai^dd cosine functions in the time and frequency domains*) 
(* BEGIN INITIALIZATION *) 



<* END XNITXAX»XV AT ION ») 

< * Program * ) 
Beg/inFackage ["Convolve * w ] 

Ko«1Cohvo1vq: :usage = "NoelConvolve [a,t»] perfonnfl a convolution, a ia a vector, b ia 
is rifltricted to length o£ a being ffreatex: ttxan length of to- e.g*. 
NoolConvolvB [ <al, a3 r a4 , a5 ) , ( {bl, cl} . {b2 , c2 ) ) ] ■ 

TeneorConvolve: :us&ge = "TanaorConvolva [a, bj performs a. convolution, a is e. vector, 
Teneorconvalve [ {al f a2 r a3 . a4 « a5} , { (bl, cl) , {b2, c2} } ] . Hee alpo NoalConvolve « 

NapConVO 1 ve I ; upagfi = 

"MapConvolve [Input . filter] performs vector convolution only. It la memoxy efficient. 
Begin [ " * Private * " J 

Noalconvolva [ f irst_„ BecondL.] x - Block [ <xl, x2 , x3 J , 
xX => Ldngtn [ second] t 
x2 = Table [Q, {i, 1, xl-1) 3 r 
X3 - Join[x2 f f irat, ac2] ; 

x4 = Partition[jK3.xl,l] ; , p 

Kap[ReverBe[#] - aecondfr, sc4j] /* *.engt:n[f irst J >= Lengthtsecond] 

HapConvolve [input_ # £iltor_] ; — 
Mbdule[{i / xl r x2 r x3), 

3cl b Table [ O , { i , 1 # Length [filter ] )]; 

Jt2 [a_3 := {xl ~ tfoin[Drop [xl, 1] „ ta>] ) ; 

ac3 = Reverse [filter] j , 

Map[(K3 . x2[»])&. Join[input,Tabla[0,{i,l,Length[filter]-l)] ] ]]* 

Ten&orconvolve [input_ r filter^] 2= Transpose t Map [HapConvolve [input , #] &. .Transpose [fil 

EruU] 

Protect [NoelConvolve , TenaorConvolve , HapConvolve] 
Null 

Enci P ackage [ ] 



11-03-98 16:41 00 44 1276 677720 

11-M/^98(WED) 16:54 NMP PATENTS UK 



P. 57 R-780 Job-222 

FAX: 00 44 1276 677720 P. 057 



bessel. m 



mwifilter s* dlordar] [0] ; 

denomfilter = Apply [Time a . tempf iltftr] ; 

roots ~ Apply[J©In,Map[S©lveC#rfl]*»Ma.p[# ===Ofc # temp«ilter] ] } ; 
ApplyCPlua. numf ilter/Dtaenomfilter, fil EMe t) /. roots] 
j / j 0 <- order < 15 

BesselFilter [order J [B&ndPass [Lower.. ,0pper_3 3 [In^uleeRoBponse] [t_] : = 
Block [ {se, fcemp£ilter,dexiomf ilter,num£iiter, roots > , 

tempfilter = (a - ^/ (d [order] [01 * ( 1/order) ) ) /. 

N[ Solve [Bessel [order] [as] ~ 0, ho)]; 
toif^filter o tempfilter/- s -> (s*2 + lower Upper ) /a / (Upper - Lower) ; 
tampfilter = simplify Co tempfilter] ; 
numfiltdr = d [order] [0] s A orderj 

denonif liter = Apply [Times, tempfiilter] t 

roots = Apply [tfoIa.MaptB©lve[#,£j]*.Kap[# B =Ot.tempfiltar] n ; 
Apply[Plua, muafilter/D [dencun£llter, s] E* (s t) /- roots] 
]/;0 <=» order < IS * 

GraupDelay [Fzreqaeacyileeponfl«„l [o_] [AngularFreauency_3 
Block[{tempfl.tampw, tl,t2.t3, t4» tampH), 

tl = FrequeneyReBponflA / . s -> temps; 

tl = tl /. temp a -> I tempw; 

tl - Together [tl] j 

t2 = Qcpand [ Denominator [ tl ] ] ; 

t2 - Distribute [Conjugate [t2] ] ; 

t2 - Map [Die tribute Times] &, t2] ; 

t2 = t2 /- Conjugate [tempx_J -> tempx; 

t2 = Expand [Numerator [tl] t2] ; 

t3 ■ DiBtributa[Rd[t2]]; 

t3 = Nap [Distribute [# r Times] & , ] ; 

t3 = t3/- RS [X_] -V x; 

t4 = ArcTan [Simplify [Expand [-1 ( t2 - t3)]]/t3); 

t4 o ~-D[t4, tmnpw]j 

t4 / - tempw — > AnguX arPr*quency 

] 

GroupDel&y: : usage a 

M GroupDelay [FrequencyHQaponfleFuac tion, Regues CedFrequency] gives the group delay 
given the f re<juency response in terms of £ at the requested frecxuency - Note the group 
delay is usually evaluated at O £or a, lorwpase filter.ft.g if 

rreouencyRafiponsofPaxPuiicjtion = 
Fuac ti on [£. Evaluate [BesselFilter [6] [LowPaee3dB[100] ] [FraQpieneyReeponfle] [2 Pi i f ] ) ] 
then GroupDelay (FrequencyRegponfleFunctimi. ££] can be ploced to show how wide the 
linear region is." 



GroupDelay [FrequencyRaflponsGFunctlon,, Requep tedFrequency^] : = 
Modulet{ £l,ael,x2,x3 > , 
xl = 

ArcTan I Camp 1 exExpand [ Im [ FrftquencyRedpons eFunc t ion [ f 1 ] ] / Ha I Fraqu oncyRespons aFunc 1 1 on [ £ 1 
1]]]* 

x2 - Evaluate [ -D [xl , £ 1] ] ; 

(x2/. £1 — > RetzuestedPrequency) / (2 Pi)]; 



(* Group Delay calculations -*> 

In[40j :=Table[N[^lB] . <W„ 2 Pi 390000,2 Pi 1010000.2 Pi 5000)] 

Out [40) « {-0,00000277055, -0.00000275661. - 0 - 00 000274284 , -0.00000272905, 

> -0. 00000271467) 

Xn [41) : ^ 1000000 2 Pi *&40 

Oufc[41)- {-5.5411 Pi, -5.51322 Pi, -5.46569 Pi, -5.4581 Pi, -5-42933 Pi) 
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1 outC42l- -0-1X17C9 Pi 
Xnt43l == * 360/12 Pil 
Out 143}= -20.1194 

*) 

~ X1 * " Block t { a. B a . teapfllter. *oot ? . l?;?'-^ '"is] 1 i 

flections - Bae tt°S? ',' % "> 2/3 (1 -">/( 1 * i U 

SiiilllilSii = ?: - 

{ acoef f i e i ent p,bcoef£ici ant d > 1 

Output [section! lFilt*rSection_, TflntntttaacniMieo ls = Blook[U«roa ( .polaa, 

zbIob = £«wtl»[Plltarsa=ti.oat£lin - 1 . 



output - singleOutput ^ixceir^^L.^, «T TTTTT - 

outiutfloaixenoa - Join [out puts e«monce, < output > 1 , 
i^put^^nce, = Drop [inptitat.ciafinco, 1) i 
Output[Bocttonl [FilterSection, 
a ixiiti*l.oixx±tiAl] [input^aquenco] [output^o^ce] 1 /- 
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BinitiaX. 1 [Inputseijutw^«-i inputsequenco .= < > 

mputS^ce Filter « (> 

*. [Drop [Pilter, 

[inputflequancollOn /• 

inputflftGwence >r a i I1 itial = ^Initial; 

^initial = Binitialj 




output soqpieTicA 1 

( + T«Sta 

0utpu u,m tBr nm^tco 8tl ooo,o 1 cT /185 * co.cxoooooo * ./ «i . o.ioiin 

aMroot8 = Ma p[^ 0 tioa[x.Pi«ritou t «[K.E«u B in.I'" t " bu " lArBl 

azrgrooca ™* 1 niBtribute [Abs [reocs , or] ] 

ab^ot* - Map [Fusion [-.niofc*""^ 1 ■ DlBt 

Plot , Aha t flit**- IX «1 1 1 . f». 0 . 5O00> 1 

*) 

EndPackagetl 



